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PREFACE 


This  project  report  is  based  on  five  lectures  given  by 
Dr.  I.  N.  Sneddon  at  North  Carolina  State  College  in  March,  1963. 
The  results  of  the  research  reported  here  have  been  successfully 
applied  to  the  solution  of  certain  crack  problems  in  the  mathe¬ 
matical  theory  of  elasticity;  some  of  these  applications  will 
appear  in  a  subsequent  project  report  and  also  are  expected  to 
appear  in  appropriate  journals. 

Copies  of  this  report  are  being  distributed  as  directed; 
this  project  is  sponsored  by  AFOSR,  ARO,  and  ONR  through  the 
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DUAL  SERIES  RELATION 

1.  Introduction, 

The  present  report  is  based  on  a  series  of  five  lectures  given  in  the  Mathematics 
Department  of  North  Carolina  State  College  in  March  1963 »  The  aim  of  theses  lectures  was 
to  present  a  connected  account  of  some  recent  researches  on  dual  series  relations,  in  par¬ 
ticular  of  some  work  done  in  the  University  of  Glasgow  by  the  author  and  Dr.  R.  P.  SrivaitaV. 

The  name  dual  integral  ecuations  seems  to  have  originated  -with  E.  Co  Titchmarsh  who 
applied  it  to  a  pair  of  equations  of  the  type 

nfXi 

j  $(K)  «(0  K(x,  0  dK  ~  f(x),  0  <  x  1  (1.1) 

c 

■  '!»(?)  K(x,  o  d?  =  g(x),  x  1  (1.2) 

**  r, 

in  which  the  weight  function  u(^),  the  kernel  K(^,x)  and.  the  free  terms  f(x),  g(x)  are 
all  prescribed  and  the  object  is  to  determine  the  unknown  function  A  survey  of 

methods  of  solution  of  equations  of  this  type  was  given  recently  (Sneddon,  1962). 

In  this  report  we  are  concerned  with  the  series  analogue  of  these  equations,  i.e.  with 
the  problem  of  determining  a  sequence  of  constants  {an‘>  satisfying  the  dual  relations 

00 

!  an  u(Xn)  K(x,  \n)  =  f  (x) ,  0  s<  X  <c  (1.3) 

OD 

Y  an  K(x,  \n)  =  g(x),  c<  xi  1  (1.4) 

X: '  I 

where  the  weight  function  w(Xn) ,  the  kernel  K(x,  X^)  and  the  free  terms  f(x),  g(x)  are 
all  prescribed,  and  iX  )  is  the  sequence  of  positive  zeros  of  a  given  transcendental 
function  j(X)»  Relations  of  this  type  seem  to  have  been  discussed  first  by  Cooke  and 
Tranter  (1959)  who  discussed  the  case  in  which  w(Xfi)  =  XnP,  K(x,  X^)  =  Jv(Xnx), 

j(X)  =  jv(x). 

The  two  types  of  dual  relation  can  be  described  within  one  framework.  Suppose  that 
the  problem  is  that  of  determining  a  function  \j>  defined  on  a  set  I  and  that 

L1  (x,0  tJ>(0  =  f X(x) ,  x  t  J1  (1.5) 

I ’2  (x,0  *  (0  =  f2(x)»  x  e  J2  (1*6) 

where  the  linear  operators  and  are  defineu  on  J  x  I  where  J  =  An 

analysis  of  this  general  problem  has  recently  been  given  by  W.  £.  Williams  (1963)  but  we 
shall  not  discuss  it  here.  Our  concern  is  with  the  solution  of  dual  relations  of  the  type 
(1.3),  (1.4)  for  special  values  of  the  weight  function  and  the  kernel. 
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The  work  described  in  this  report  arose  out  of  the  analysis  of  certain  physical 
problems  so  we  begin  in  (.  §  2)  by  reviewing  some  of  the  situations  in  which  dual  series 
relations  occur  in  the  analysis  of  mixed  boundary  value  problems  in  mathematical  physics. 
The  solution  of  these  dual  series  relations  involves  complicated  mathematical  analysis. 
With  a  view  to  shortening  the  proofs  by  separating  out  some  frequently  occurring  pieces 
of  formal  calculation  we  separate  out  (in  f3)  some  of  the  basic  mathematical  techniques. 
In§?4,  5  we  consider  the  various  kinds  of  dual  series  relations  in  which  the  kernel 
K(x,  X  )  is  a  Bessel  function  of  the  first  kind.  The  relations  are  essentially  of  two 
types:  in  the  first  (§4)  the  function  j (X)  is  the  Bessel  function  J  (X)  and  the  series 
involved  are  therefore  Fourier-Bessel  series,  while  in  the  second  (^5)  the  function  j (X) 

i 

is  the  linear  combination  XJ  (X)  +  HJ  (X)  and  the  series  involved  are  Din!  series.  In 

y  V 

§56,7  we  discuss  dual  relations  involving  trigonometrical  series  which  are  of  use  in 
the  solution  of  certain  two-dimensional  boundary  value  problems.  In  the  next  section 
(§  8)  we  discuss  a  class  of  dual  series  relations  involving  series  of  associated  Legendre 
functions,  using  the  method  of  Collins  (1961) .  Series  relations  of  this  kind  which  arise 
in  the  analysis  of  problems  relating  to  spherical  caps  form  a  particular  case  of  the 
problem  of  solving  certain  dual  series  relations  involving  series  of  Jacobi  polynomials. 
We  conclude  our  treatment  by  outlining  (in  48)  Srivastav’s  method  of  solving  equations  of 
this  type. 


...  The  Occurrence  of  Dual  Series  Relations  in  Mathe’natical  Pir/uicu. 

We  begin  by  considering  some  of  the  circumstances  in  which  dual  series  relations 
arise  in  the  analysis  of  boundary  value  problems  in  mathematical  physics. 

2.1 .  Electrified  Disk  in  a  Grounded  Cylinder. 

Suppose  that  an  electrified  disk  is  situated 
in  a  plane  normal  to  the  axis  of  an  infinite 
circular  cylinder  its  centre  lying  on  that  axis. 

The  cylinder  is  assumed  to  be  grounded,  and  the 
potential  of  the  disk  is  prescribed.  We  take  . . 
the  origin  of  coordinates  at  the  centre  of  the  disk 
and  the  z-axis  along  the  axis  of  the  cylinder.  For 
convenience  we  can  take  the  unit  of  length  in  our  problem 
to  be  the  radius  of  the  disk;  we  may  then  take  the 
radius  of  the  cylinder  to  be  a  >  1 .  The  problem  is  to 

find  the  potential  V  of  the  electrostatic  field  in  the  p’-!?-  1 

interior  of  the  cylinder.  If  we  introduce  cylindrical 

coordinates  (  p  ,  p  ,  z)  then  V(  p  ,  p  ,  z)  must  be  a  solution  of  Laplace's  equation 


3  Y  ^  1  37  ^  1  9  V  ^  3ZT 

dp 

where  V(  p  ,  p  ,  z) 


3 p 2  p3p  p2  dp2  3zZ 


(2.1) 


►  0  as  |  z  |  fc  oo  and 

V(p,  p  ,  0)  =  f  (  p  ,  p),  0  $  p  i  1, 

and  V(a,  p  ,  z)  =  0,  |  z|  >  0.  (2.2) 

Using  the  superposition  principle  we  see  that  if  we  can  solve  the  boundary  value 
in  which 

v(  P  >  <t>  ,  0)  =  f„(p)  cos(v  p  +  a  v),  0sp<1,  (2.3) 

we  can  solve  the  boundary  value  problem  in  which  V(p,  p  ,  0)  =  f  (  p  ,  p  ) 

Instead  of  solving  the  Dirichlet  problem  for  the  space  between  the  disk  and  the 
circular  cylinder  we  can  reduce  the  problem  to  a  mixed  boundary  value  problem  for  the 
semi -infinite  cylinder  z.*0,  0«p?a  when  the  boundary  conditions  (2.2),  (2.3)  are 
supplemented  by  the  condition 

97  =  °»  z  =  °»  1  <  P  <  a.  (2.4) 

The  function 

V(  P  ,  p,  z)  =  cos (vp  +  a„)^  X”1  an« 

n*i 


-X  z 


(2.5) 
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satisfies  the  harmonic  equation  (2.1 );  it  satisfies  the  condition  V  — 0  0  as  z  CD 
and.  the  boundary  condition  (2.2)  provided  that  X  ,  Xg,  ^  ,  ...  are  chosen  to  be  the 
positive  roots  of  the  transcendental  equation 

J,(Xa)  =  0.  (2.6) 

To  satisfy  the  boundary  conditions  (2.3),  (2.4)  we  must  choose  the  constants  ,  a^, 
a3 ,  ...  to  satisfy  the  relations 

OO 

^  ^  n  a  ^  (p  )  ,  0^p<1  (2.7) 

n  =  i 


OO 


W)  -  °- 


1  <  p  «  a. 


(2.8) 


In  particular  if  the  prescribed  potential  of  the  disk  is  symmetrical  about  the 
z-axis,  i.e.  if  V(  p  ,  6  ,  0)  =  f (p  )  ,  for  0  £  p  <  1 ,  we  may  take  v  =  0  in  these 
equations  to  obtain  the  pair  of  equations 


oo 

Y, 


0  «  p  <  1  , 


(2.9) 


00 


1  <  p  «  a. 


(2.10) 


1 1  ~  1 

where  now  the  f\  p J  are  the  positive  zeros  of  the  function  J^(Xa)t 

Pairs  of  equations  of  the  type  (2.7),  (2.8)  or  the  type  (2.9),  (2.1 0)  in  which 
the  sequence  of  constants  JW  is  determined  by  a  pair  of  equations,  one  of  which 
is  valid  in  one  segment  of  the  positive  real  axis  and  the  other  of  which  is  valid  over 
another  segment,  are  called  dual  series  equations. 


2,2.  The  Beissner-Sagoci  Problem  for  a  Semi -Infinite  Cylinder. 

The  standard  Reissner-Sagoci  problem  is  that  of  determining  the  components  of 
stress  and  displacement  in  the  interior  of  the  semi-infinite  elastic  solid  z  >  0  when 
a  circular  area  of  the  boundary  surface  z  =  0  is  forced  to  rotate  through  an  angle  a 
about  an  axis  which  is  normal  to  the  undeformed  plane  surface  of  the  solid;  it  is 
assumed  that  the  part  of  the  boundary  surface  which  lies  outside  this  circle  is  free 
from  stress.  (Of.  Sneddon,  1951*  p.500).  Reoently  Sneddon  and  Srivastav  (1963) 
have  generalized  this  problem  to  the  case  in  which  the  semi -infinite  elastic  body  has 
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a  cylindrical  boundary  whose  cross  section  is  a  circle  of  radius  a. 

It  has  been  shown  by  Reissner  ( 1 937)  that  in  this  case  the  only  non-vanishing 

component  u  ,  -  where  we  use  cylindrical  coordinates  (  p  ,  9 ,  z) ;  all  the  components 
9 

of  the  stress  tensor  vanish  except  a  .  and  a  ,  which  are  related  to  u  .  through 

e  z  <p  p  <p  9 

the  equations 


3u. 

as9  =#i  Tz  > 


V  =  PP 


9j  (jt\ 

3p  l  P  /  ' 


(2.11) 


where  p  denotes  the  rigidity  modulus.  To  satisfy  the  equation  of  equilibrium 

do  .  3a,  2a. 

Bp  +  dz  +  p  ~  u 

we  must  choose  a  form  for  u  ,  satisfying  the  equation 

9 


1  9  u9 

)2  p  dp  p2  dz 2 


(2.12) 


The  boundary  conditions  of  the  problem  are  simplified  if  we  choose  the  radius  of 
the  prescribed  circle  to  which  the  rotation  is  applied  to  be  our  unit  of  length.  On 
the  plane  boundary  z  =  0  we  then  have  the  boundary  conditions 


u,  =  f(p  ),  0  S  p  <  1  ; 


=0,  1  <  p  <  a. 


(2.13) 


v= 0 


?iR.  2 


where  the  function  f (p  )  is  prescribed.  In  the  case 
of  most  immediate  physical  interest  f ( p )  -  a  p  , 
where  a  is  a  constant,  but  it  is  of  value  to  derive 
the  solution  in  the  general  case. 

In  the  first  instance  we  assume  that  the 
cylindrical  wall  of  the  cylinder  is  rigidly  clamped; 
i.e.  we  assume  that 


u9  =  °»  P  =  a>  z  >  °* 


(2.14) 


We  also  assume  that  u  ^ ,  a  z<t>’ap<t>  ^en<*  zero 
as  z  oo  .  (Of.  Fig.  2). 


If  we  observe  that 


(JL  +  12.1. 

\3p2  P  3P  P2: 


=-2. 


dz  J  dp  dp 


we  see  that 


9  3p 


(2.15) 


is  a  solution  of  equation  (2,12)  provided  that  9 
is  a  harmonic  function.  In  this  case  the  expressions  for  the  non-vanishing  corapon- 
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enta  of  the  stress  tensor  assume  the  forms 


azrp  p  dpdz  ’ 


3P*  =  12 p 


jl(  1 ) 

dp\~r)p) 


If  we  7/rite 


ifi(p  ,  z)  =  -  /  X  a  e  n  J  (X  p) 
•  ?  /  «  n  n  or. 


(2,16) 


(2.17) 


in  equations  (2.* 5)  and  (2.1 6),  where  X  ,  Xg,  ...,Xn,  are  the  positive  zeros  of 
J  (\a) ,  then 

1  _X  z 


u*  L 


X  ae  n  J  ( X  p) 

n  n  in 


(2.18; 


C  =  ~  u 
z 


-  X  z 

ae  j(Xp), 


(2.19) 


so  that  equation  (2.14)  is  satisfied  and  the  mixed  boundary  conditions  (2.1 3)  will  be 


satisfied  if  we  choose  the  constants  &n  to  satisfy  the  dual  series  equations 

OO 

X';  an  J,(Xnp)  =  f(p),  0«p  <  1, 

n  =  1 

00 

£  .,/(>,,)  =  C,  1  <P  «  ». 


(2.20) 


(2.21) 


sre,  it  7/ill  be  remembered,  the  w  are  the  positive  zeros  of  J^(Xa). 

If,  instead  of  the  boundary  condition  ( 2 . 1 4)  ,  we  assume  that  the  cylindrical 


boundary  is  stress-free,  i.e.  if  we  assume  that 


'tT* 


%4>  =0>  p  =  a>  z  >  °* 


(2.22) 


Oy 

7  r 


J 


it  is  necessary  to  consider,  not  a  semi -infinite 
•  CT,  =  0  cylinder,  but  a  very  long  one  of  length  6  >>  a. 

Pt 

We  assume  that  the  base  of  this  cylinder  has  zero 

c  displacement,  i.e.  that 
2*0) 

-  £ y  =  0,  Z  =  6  ,  Os;  p  =;  a.  (2.23) 

z'  The  conditions  (2.13),  (2.22),  (2.23)  can  be  realized 
/  by  considering  the  distortion  of  a  cylinder  which  is 

rigidly  attached  to  a  fixed  rigid  foundation  z  =  6 
and  is  deformed  by  the  application  of  a  torque  to  a 
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circular  patch  of  its  other  plane  surface  (s.g.  through  the  rotation  of  a  rigid  disk 
attached  to  it)  while  the  remainder  of  its  surface  is  free  from  stress  (Cf.  Fig.  3) 
If  we  assume  the  form 

sinh  X  ( 6  -  z) 
n 

sinh  X  6 
n 

n  =  i 


for  the  displacement,  then  we  see  from  equations  (2.1 1 )  that  the  stress  components 
are  given  by  the  equations 


EcoshX  (  6  -  z) 

an  sinh  X  n&  J,  (  V  }  ' 


°  A.  =  ~  V  a 

z<p  0 


mm  sinh  X  (  6  -  z) 

O*  =  -  )  a  - t- r? — ? -  J  (  X  p) 

/  (  n  sinh  X  6  2  n 

n  =  i 


(2.25) 

(2.26) 


It  fellows  immediately  from  these  equations  that  the  boundary  condition  (2.23)  is 

satisfied  and  that  the  boundary  condition  (.2.22)  will  be  satisfied  provided  the  X  's 

n 

are  the  positive  zeros  of  Jz(Xa).  Also,  since  on  z  =0, 


a'nJ,(XnP) 
n  =  i 

oo 

a z<}>  =  -paoP  “  ^  ^  an  coth(Xn6)  J  (Xnp) 

n=  i 

it  follows  that  the  boundary  conditions  (2.1 3)  will  be  satisfied  if  the  coefficients 
a  satisfy  the  dual  series  relations 


OO 

sp6  +  ^  X"1  a n  J,(Xnp)  =  f(p),  0  s  p  <  1 


(2.27) 


aQp  +  a^  coth(Xn6)  (X  np)  =  0,  1<p  s  a. 


The  last  pair  of  equations  is  the  pair  appropriate  to  the  deformation  of  a 
cylinder  of  finite  length  6  .  If  we  are  considering  a  very  long  cylinder  (5  »  a) 

then  we  may  replace  the  factors  coth(Xn§)  occurring  in  the  series  on  the  left-hand 
side  of  equation  (2.28)  by  unity.  For  a  long  cylinder  the  problem  is  therefore  re¬ 
duced  to  the  solution  of  the  dual  series  relations 
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a  p6  + 
0 


a0  + 


^  X"n  an  J1  (Xnp)  =  fk)»  0  «  P  <  1  » 
r»*i 

oo 


(2.29) 


(2.30) 


where  8  is  a  known  constant  and  the  X's  are  the  positive  zeros  of  J  (Xa). 


2.3 »  Problems  in  the  Conduction  of  Heat. 


Si  milar  equations  arise  in  problems  concerned  with  the  conduction  of  heat  in 
cylinders.  If  6  (  p  ,  <p  ,  z)  la  the  deviation  of  the  temperature  at  a  point  with 
cylindrical  coordinates  (  p  ,  <f>  ,  z)  from  a  standard  temperature  0Q  then  it  is  well- 
known  that  in  the  steady  state  0  must  be  a  harmonic  function  in  the  region  considered, 
so  that  a  solution  appropriate  to  the  semi -infinite  cylinder  0  ?  p  s  a,  z  >  0,  in 
which  there  is  an  axisymmetric  flow  heat  is 


-X  z 

XnanJ„(pXn)e  " 

n  n  o  n 

n  =  i 


(2.31) 


From  this  equation  it  follows  that  the  rate  of  flow  of  heat  per  unit  area  across  the 
plane  z  =  0  out  of  the  solid  is 


(*“*),.„  *  -k  £  %  W- 


(2.32) 


where  k  is  the  conductivity,  while  that  across  the  cylindrical  surface  p  =  a  out  of 


the  solid  is 


('“^La  Z  a.e’Vj>X") 


(2.33) 


The  constants  Xn  are,  as  yet,  unspecified.  Their  values  will  be  determined  by 
the  boundary  conditions  on  the  surface  p  =  a.  There  are  three  important  cases  to 
be  considered. 

(i)  In  the  first  case  we  assume  that  the  temperature  of  the  surface  p  =  a  is 
maintained  at  the  constant  value  0Q  so  that  0(a,  z)  =  0.  From  equation  (2,31 )  it 
follows  that  in  this  case  the  XJs  are  chosen  to  be  the  positive  zeros  of  J  (Xa)„ 

If  the  plane  surface  z  =  0  is  heated  in  such  a  way  that  the  temperature  is 
prescribed  over  the  circle  0  <  p  <  1  and  the  remainder  of  the  surface  is  insulated  then 


9 


we  obtain  the  pair  of  dual  .integral  equations 


OO 

j  \~’f  a  J  (p  X  )  =  f (p),  0$  p  <  1 

/  n  n  o  n  r '  * 

n*  i 

OO 

i<p<a> 


(2.34) 


(2.35) 


where  the  function  f (p)  is  prescribed. 

On  the  other  hand  if  the  surface  is  heated  in  such  a  way  that  there  is  a 
prescribed  flux  of  heat  into  the  circle  0  <  p  <  i  and  the  remainder  of  the  surface 
is  maintained  at  the  standard  temperature  d then  the  a  must  satisfy  the 
equations 


OO 

Za  J  (pX  )  =  g(p),  0«p  <  1 

n  c  n 

n=  i 

OO 

LX"1  a  J  (pX  )  =  0,  1  <  p  £  a. 

n  n  on'  '  r  * 


(2.36) 


(2.37) 


where  the  function  g(p)  is  prescribed. 

Equations  of  a  more  complicated  type  arise  if  we  have  radiation  from  the  surface 

z  =  0  into  a  medium  maintained  at  a  constant  temperature  Q  . 

0 

By  Newton's  law  of  cooling  we  then  have  that  on  z  =  0 


W. 


H6(p  ,  0) 


where  H  is  a  constant.  From  equations  (2.31 )  and  (2.32)  we  then  find  that 

OO 

(h«  -  k  |f)  =*  Yj  a"(l  *  h’7')J.('>\.) 

where  h  =  B/k. 

Hence  if  the  temperature  is  prescribed  over  the  circle  0  <  p  <  1  of  the  surface 
z  =  0  and  over  the  rest  of  the  surface  there  is  radiation  into  a  medium  maintained  at 
temperature  Qq  ,  the  constants  an  occurring  in  the  solution  (2,31 )  must  satisfy  the 
dual  series  equations 


00 

Yj  x"n  an  J0(v)  =  0«P<1, 


(2.38) 
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oo 

l 

n- 1 


a  (l  +  hx’1  )j  (x  p)  =  0,  1  <  M 

n  n  a  n  ' 


(--39) 


On  the  other  hand  if  the  heat  flux  is  prescribed  over  the  circle  z  =  0,  0  ?  p  <  1  and 
we  have  the  radiation  condition  over  the  zone  1  <  p  $  a,  we  see  that  the  equations 
determining  the  coefficients  a  are 


OO 

£ 


n«  i 


oo 

£ 

n~i 


an  =  s(p),  0.  p  <  1 

(l  +  h Xn ’)  a n  Jc(Xnp)  =0,  1  <  P  «  a 


(2.40) 

(2.4 1 ) 


(ii)  In  the  second  case  we  assume  that  the  cylindrical  surface  p  =  a  ±3 
insulated  against  the  flow  of  heat.  Prom  equation  (2.33)  it  follows  that  the  \  ?s 
in  thus  case  are  the  positive  zeros  of  J  (Xa).  We  again  get  boundary  value  problems 
of  the  types  considered  in  (i) y  The  dual  series  relations  are  exactly  the  same  as 
those  listed  above  except  that  the  Xn  are  now  the  positive  zeros  of  J  (  Xa) . 

(lit)  Finally  if  there  is  radiation  from  the  cylindrical  surface  into  a  medium 
maintained  at  the  fixed  temperature  6  we  .have,  as  a  consequence  of  Newton's  law  of 
cooling,  the  equation 


(«  ♦*§£)  -0- 

'  s  p  =:  a 

Substituting  from  equations  (2,31 )  and  (2 .33)  we  find  that 
/  00 

*kf f)  =  *  Y  ■►v.’tv?}' 

^  P  —  3 


(2.42) 


-X.  z 
n 


from  which  it  follows  that  the  conditions  (2 *-4-2)  is  satisfied  if  we  choose  the  X 
occurring  in  the  solution  (2»3l)  to  be  the  positive  roots  of  the  transcendental 
equation 

hJQ(Xa)  +  X,F(Xa)  =  0.  (2.43) 

Dual  series  involving  trigonometric  series  arise  in  the  solution  of  problems 
concerning  the  conduction  of  heat  in  the  long  strip  0  $  x  s_  n  ,  0  ^  y  5  4  If  we 
suppose  that  the  temperature  deviation  d  (x,  y)  is  prescribed  over  the  segment 
0«x<c,y  =  0  and  that  the  remainder  of  the  boundary  of  the  strip  is  insulated 

the  flow  of  heat  across  it  then  at  is  readily  shown  that  the  solution,  of  the 
problem  will  be 
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e(x,  y)  =  ?a0(6  -y)  +  ^  n'1ancos(nx)-22^^~:-^,  (0  <  x<  7r,0<y?5),  (2,44) 

n  =  i 

provided  that  we  can  determine  a  sequence  of  constants 
series  relations 


oo 


i*„6 

n  ’a  cos(nx)  =  f(x),  0  $  x  <  c, 

» 

OJ 

v — / 

n  =  i 

oo 

?  ao +  y 

a  tanh(n6)  cos(nx)  =0,  c  <  x  5  jr, 

(2.46) 

n  =  i 


where  the  function  f(x)  is  prescribed.  If  6  >>w  we  can  replace  tanh(n6)  on  the 
left-hand  side  of  equation  by  unity  to  obtain  the  dual  series  relations 


M 


to  satisfy  the  dual 


1 

2 


a  6 
o 


i  a 


co 

^  n_1a  ^  cos(nx)  =  f(x), 
n  *  i 

OO 

V  a^  cos(nx)  =  0, 


n  =  i 


0  ?  x  <  c. 


c  <  x  s  w. 


(2.47) 

(2.W) 


On  the  other  hand  if  the  conditions  on  the  sides  x  =  0,  x  =  v  are  replaced  by 
the  conditions  A  (0,  y)  =  A  (  tt  ,  y)  =  0  the  solution  appropriate  to  the  semi-infinite 
strip  0  «  x  <7r,y*0is 


OP 

A(x,  y)  =  ^  n"1  an  sin(nx)e "nu  , 


(2.49) 


n=  i 


N 


where  the  constants  a  satisfy  the  dual  series  equations 


^n-  a  n  sin(nx)  =  f  (x)  ,  0  *  x  < 

n»i 

oo 

y  an  sin(nx)  =  0, 


(2.50) 


n*i 


c  <  x  $  no 


(2.51) 
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2.4-.  Some  Boundary  Value  Problems  for  a  Spherical  Cap. 

Dual  series  relations  involving  F'ourier-Legendre  series  arise  in  the  analysis 
of  some  electrostatic  and  hydrodynamic  boundary  value  problems  for  a  spherical  cap. 

The  first  problem  vre  consider  is  that  of  determining  the  electrostatic 
potential  u  due  to  a  thin  spherical  cap  maintained  at  a  prescribed  potential.  If  we 
use  spherical  polar  coordinates  (r,  6  ,  <p)  referred  to  the  centre  of  the  sphere  as 
origin  and  the  axis  of  symmetry  of  the  cap  as  polar  axis  (Oz  in  Fig.  4)  ,  we  can 

describe  the  cap  by  the  equations  r  =  a, 

0  ^  9  $  a  .  On  the  assumption  that  the 
prescribed  potential  is  symmetrical  about  the  axis 
of  the  cap,  i.e.  that  u  =  tie)  on  the  spherical 
cap,  vre  may  take  u  to  be  an  axisymmetric  solution 
of  Laplace's  equation 


Fig*  h 


02 
0  U 


0  r 


2  0u 

+  —  —  + 

r  9r 


sin 


Odd  \ 


sinfl 


3u  \ 

d  e) 


satisfying  the  conditions 

u  =  f(fl)  on  r  =  a,  0  $  9  $  a  ;  and  u  is  continuous  across  r  =  a; 

3  v/d  r  is  continuous  across  the  spherical  region  r  =  a,  a  <  8  <  n ; 
u  =  0(r”1)  for  large  r. 

A  solution  of  Laplace's  equation  satisfying  the  continuity  condition  on  u  across 
r  =  a  and  the  condition  (iii)  is 


(i) 

(ii) 

(iii) 


OO 

£a„\f)  Pn(cOSe)»  °"r  <-  a 


u(r, 6)  = 


n  =  o 


1  Pn(COS0)» 


r  <  a 


where  to  satisfy  the  conditions  (i)  and  (ii)  the  constants  a  must  be  chosen  to  be 

n 

solutions  of  the  dual  series  relations 


Lj 

n*o 


Pn(cos  e )  =  f(e). 

(0  5  9  <  a  ), 

(2.52) 

P  (cos  9 )  =  0, 
n  9 

(a  <  6  s  2jt). 

(2.53) 

t 
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A  similar  boundary  value  problem  involves  the  determination  of  the  velocity 
potential  ft  due  to  the  motion  of  a  spherical  oap  along  the  direction  of  its  axis  with 
constant  velocity  U  in  a  perfect  fluid  at  rest  at  large  distances  from  the  cap.  If 
we  use  the  above  coordinate  system  then  t  =  §  (r,  $  )  must  be  an  axisymmetrical 
solution  of  Laplace's  equation  such  that,  for  large  values  of  r,  ft  -  o(r"2)  and 
rd$/ 9r  is  continuous  across  the  sphere  r  =  aj  further  it  must  satisfy  the 
condition 

r|*:  =  -UaooBfl,  r  =  a,  0  <  ft  *  0  (2.54) 


while  on  the  spherical  region  r  =  a,  a  <  ft  *  w  ,  ft  moat  be  continuous  (Lamb,  1955, 
p.l60).  An  s-rlsj  man  trie  solution  which  is  0(r-1)  for  large  r  and  has  r  S  i/3  r 
continuous  across  r  =  a  is  given  by  ■Hi«i  equations 


*(r,e)  = 


2^(n*l)a#^J  P^  (  cos  ft  )  ,  (r  <  a), 
n-i 

’!/*"(*)  (r  >  a), 

n*  t 


The  condition  (2.54)  and  the  continuity  of  ft  across  r  *  a,  a  *  0  *  r  are  satisfied 
if  the  constant*  a^  satisfy  the  dial  series  equations 


l)a  P  (ooa  ft)  s  -  Ua  cos  ft,  0  <  ft 

n  n  w 


<  o. 


a»i 

m 

Z(2n  ♦  l)an  pn(°°*  O  *  0»  o  <  ft  <  w. 

n»i 


(2.55) 

(2o56) 


On  the  other  hand  if  me  suppose  that  the  oap  aovea  through  the  fluid  in  the  direction 
4>  =  0  perpendicular  to  its  axis  of  symmetry  with  uniform  velocity  U.  The  velocity 
potential  ft  (r,  ft  ,  $)  for  the  motion  must  satisfy  Laplace's  equation 

r*  ill  .  a,  81  .  _1_  A.  ii) 

3r  dr  sin  ft  3  ft  ^  3  »J  *in*ft  30* 

be  such  that  ft  =  0<r_I)  for  large  r,  ft  and  3  ft/3 r  are  continuous  over  the  surface 
r  s  a,  a  <  ft  <  2f  while 


*•  =  -  U  a  sin  ft  ooa  . 


(2.57) 
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It  is  easily  shown  that 


|  y^(n  +  2)a^  (cos6)cos<£,  0  «  r  < 

*(r  ,e,4>)  =  J 


Z(n  -f-  1 )  a  /  a  \  T  (cos  6  )cos  <f> 

n  f  -  J 
n=o  '  ' 


where  the  coefficients  a  satisfy  the  dual  series  relations 

n 

oo 

yi  (n  +  1  )(n  +  2)  (cos  0  )  =  -Ua  sin  Q  ,  0  <  0  <  a 

n*i 

oo 

£(2=0)..  Tn^  (cos  0  )  =  0,  a  <  0  <  7T  . 

n  =  i 


3.  Mathematical  Preliminaries . 

The  solution  of  the  various  types  of  dual  series  relations  which  we  have  derived 
in  §2  is  not  a  simple  matter.  The  analysis  is  complicated  even  in  the  simplest  case 
and  it  is  preferable  to  shorten  the  proofs  by  separating  out  various  bits  of  formal 
analysis.  This  is  what  we  shall  do  in  this  section.  We  begin  by  listing  some 
properties  of  operators  of  fractional  integration. 


The  Erdelyi-Kober 


srators. 


We  shall  make  use  of  the  Erdelyi-Xober  operators  I  as  modified  by  Sneddon 

(1962).  If  T)  and  a  are  real,  a  *  0  and  7?  >  -  £  we  define  the  operator  I  a 
by  the  equation 


11  ]- 


-2  a  -2  r\  /•* 


x  2  77  +  1  /  2  2\a  -  1  „/  \ 

u  '  (x  -  u  ;  f(u;< 

*  a 


(3.1) 


while  if  a  <  0  and  n  is  the  least  positive  integer  such  that  a  +  n  >  0  we  define 

it  by  the  equation 


4 ,  .[*»>«  -  ]  =  *-2’ ' ■ 20  - 1  a;  *2"  ♦  ■ 20  +  2n  +\.+ £o, ,] 

where  35  denotes  the  operator  defined  by  the  equation 

X 


(3.2) 
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4f  4s(x)- 


(3.3) 


We  shall  in  fact  only  encounter  the  case  in  which  -«1  <  a  <  0  when  equation  (3.2) 


reduces  to  the  simple  form 


A  much  used  property  of  these  operators  is  that 

f(u)*  X]  =Cx2^Ir,  +  (9,a[f(u)’*  3 


(3.4) 


We  shall  also  require  the  fact  that  X  has  an  inverse  defined  by  the  equation 

Tl  t  a 


and  the  property 


I"1  =  I  , 

r?  ,  a  T)+  ay  -  a 


X  .y  I  «  —  X  .3  • 

n  ,  a  Tl  +a,p  T)  ,  a  +P 


(3.6) 


(3.7) 


We  also  note  that  1^  q  is  the  identity  operator. 


When  the  function  to  which  the  operator  1^  a  is  applied  is  a  function  of  mere 
than  one  variable  we  shall  modify  the  notation  of  (3.l).  For  example  we  shall  use 


,  a  ff(u,  t);  u  xj  to  mean 

~T(a)  -P*2"  ■"1(x2-u!)“-1  f(»,  t)dn. 


It  considerably  simplifies  the  writing  down  of  our  results  if  we  define  a 
function  K  g  (  p  ,  t;  a)  to  be  the  integral 

vt  a  »  p  1  y 


r  K,(ay) 

ip  +y  ^ 

(when  the  integral  converges).  When  it  is  clear  from  the  context  what  v  and  a  are  we 

shall  write  K  „  (p  ,  t)  for  K  .  (  p,  t:  a).  The  integral  will  satisfy  the 

correct  convergence  conditions  'at  infinity*  if  2a  >  P  +  tj  since,  in  our 

applications  p  <  a,  t  <  a  this  condition  is  automatically  satisfied.  To  ensure 

convergence  at  the  lower  limit  of  integration  we  must  have  -2  v  +  a  4.  p  +1  +  y  >  -1, 

i.e.  we  must  have  v  <  l+j(a+P+y)» 

The  following  relations  satisfied  by  K  R  (p  ,  t)  can  easily  be  deduced  from 

a  t  P  t  Y 

the  recurrence  relations  satisfied  by  modified  Bessel  functions  of  the  first  kind 


(Watson,  (1944),  p.79)  and  integrals  given  by  formulae  (6),  (7)  on  p.365  of 
Vol.II  of  Erdelyi  (l954): 
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(5’ =  1, /S,  (3.8) 

y  + 1  Co.‘).  (3.9) 

v^’ t)!  y°  *}  *  *><  <3.io) 

t+^y-J*'  t).  (3.11) 

We  shall  also  have  occasion  to  consider  the  function  defined  by 

k*  /  s  ry*l( y)  +  HS(y)  „  , 

>',Hl^y,6  lu»v)  =1  — —  —  ^(uy)l  (vy)y1+6  (3.12) 

Jo  yl’(y)  +  Hiv(y)  H  Y 


(0<u<l,0<v<l)  whenever  the  integral  converges.  Using  the  recurrence  re¬ 
lations  satisfied  by  modified  Bessel  functions  of  the  first  kind  and  the  results  i 1 
Erdelyi  (l954)  quoted  above  we  find  that 


“■if"**  X“,H,iS,y,s(u>  v)}=  ^  Kv,n,  f)  -  i,y,6.,(".  »>.  (3.13) 

*tr  JS»,S,f,y,  S(u>  v)}*TyKt,H,^,y-1(  S  *1  (u*'r)>  (3.14) 

hf,  «f *.»,/».  V,  6 <u-  -  2“’‘'“<jH,a+/J,y,6-aC*>''>  (3.15) 

V)!  T^y]  “2  5,iKv’,H^,y*a,6-i“.y)(3.l6) 


3.2.  Infinite  Series  involving  Bessel  PunotTona- 

In  the  subsequent  analysis  we  shall  encounter  functions  of  the  type  S 
(p,  t;  a)  defined  as  the  sum  of  the  infinite  series 

_  /  .  2  J  (p  X  )  Jfl(tx  ) 

Sy,a  S  y('’>tS  * )»—  )  *  "  P  "V, 

■M  n 


(3.17) 


whenever  the  infinite  series  converges .  The  infinite 


M 


sequence  U  is  formed  of  the 


positive  zeros  of  the  function  J,,  _,(aX).  In  equation  (3.17),  P  and  t  are  real 

^  "I*  I 

variables  satisfying  the  inequalities  0  $  p  ,  t  5  1  <  a  and  -(  a  +  +  2)  <  y  <  0. 

Series  of  this  type  are  considered  in  Sneddon  and  Srivastav  (l  9 63)  but  special 
cases  were  considered  earlier  by  Tranter  (l  959)  •  Tranter  defines  two  functions  by 


the  series 


/  ,  %  °  v  +2m  +  k'  s' 

k*')  =  L  7r?-77ir 

.  s  y  +  1  s 


(Xa)  Jv(“aP) 


v  >  -1 


(3.18) 


*  J 
n,  k)  =  V  - 

2  L-i  X 


jc(^a)Jv  2n  +  k^s^  ,  r ,n\ 

■5 — 5 - 1 -  ,  v  -m-n-k  (3.19) 


X  2  J2  .(aX  ) 

S_  V  +  1  s 


where  m  and  n  are  zero  or  positive  integers,  k  >  0,  and  a  >  -  1 .  It  is  easily  seen 
that,  in  the  notation  of  equation  (3*1 7)  » 


*,(  V  ,  m,  k,  p  )  =  £  a  Sv>  Vf  v  +  2m  ♦  k,  -k  ^p*  1j  a^»  (3.1 8a) 


.2(.,.,n,k)  =  i>2  Sy>  „  „2„1.k,»»2n  *  k,-2(l  >' >a)-  (3-1S*> 


We  shall  also  find  that  we  have  to  deal  with  series  of  the  type 


,♦  (u  T)  -2y  VuX.)Jy(rt.A 
"’E’^y’s  k 

)<u<1,0<T<1,  .here,  no.,  f  Xa  j  is  the  ij 
positive  roots  of  the  transcendental  equation 


(3.20) 


0  <  u<  1,0<  v<  1,  where,  now,  jX^  |  is  the  infinite  sequence  formed  by  the 
positive  roots  of  the  transcendental  equation 

X  JJ(X)  +HJV(X)  =  0,  (3.21) 

H  and  v  being  real  constants,  v  >  -  Series  of  this  type  have  been  discussed  by 
Srivastav  (1963a). 

To  obtain  an  integral  representation  of  the  series  on  the  right-hand  side  of 
equation  (3.17)  we  consider  the  contour  integral 


/  P(z)dz 

Ja 


where  the  function  P(z)  is  defined  by  the  equation 
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(a  z)  +  i  Xu(a  z) 

*■)  *=- — rtny — ^V*'1 


(3.22) 


and  the  contour  0  consists  of:- 

(^i)  the  portions  of  the  positive  real 

axis  joining  the  points 


s  +  1 


4i  It 


6  ,  X1  -  81  ;xs  +5S,  xs  +1"  6 

(s  =  1 ,2,  ...  ,  p  -  l)j  Xp  +6p,  Hj  where 

the  5  '  s  are  small  and  p  and  R  are  large 

and  such  that  X  <  R  <  X  . ; 

p  p  +  1  * 

(ii)  a  series  of  small  semi -circles 
y  (s  =  1,  2,  ...,  p)  with  equations 

|z  -  X  |  -  6  j 
1  s  s'  s' 

(iii)  a  large  circular  quadrant  |z|  =  R, 

0  *  arg  z  <  3-  jr  J 

(iv)  the  positive  imaginary  axis  from  z  =  Re  2  xv  to  z  =  6e2 

(v)  a  small  circular  quadrant  |z|  =  6,  0  $  arg  z  s  ^ir. 

In  general  the  point  z  =  0  is  a  branch  point  of  F(z).  Each  branch  is,  however,  an 
analytic  function  and  therefore  this  presents  no  difficulties;  we  choose  that  branch 
for  which  Re  |p(z)J  =  J a(  p  z)  J^(  t  z)z  *  y  for  real  values  of  z.  Using  standard 
procedures  in  the  calculus  of  residues  we  can  easily  show  that 

Sv,a,0,y(p»  t;a)  =/  Jc/px)J  (tx)xy+1dx  +  |sinji(a+/9  +  y-2v)J  y(P»  t|  «0  (3.23) 

where  the  function  K  R  ( p  ,  t;  a)  is  the  integral  defined  above. 

v>  a  t  P  »  Y 

In  a  similar  way  we  can  obtain  an  integral  representation  of  the  series 
*  ♦ 

S  „  (u.  v)  in  terms  of  the  integral  K  (u,  v)  by  considering 

v,H,/9,y,8  *  v,H,  /?,y,6 

the  contour  integral 

z)J^(u  z  )  Jy(vz)z  8  +  1  <3z 

Where  /  *  j;(z)  +  ii;(z)  +  h  j  (z)  +  it  (z) 

M  = — - - - - - - — 

2j;(z)  +  H  Jy(z) 

and  0  is  the  same  contour  as  before,  except  that  now  the  X  's  are  the  positive  roots  of 
equation  (3.2l).  We  find  that 
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St,H,/J,y,6(u'  v)  =  j/f 9(ux)  Jy(vx)  ^ 

+  ^aln^6  +  0  +  y  -  2i/)Jk*^h^  y  ^  &(u,  v)  (3.24) 

We  see  that  the  expressions  for  S  and  S*  involve  the  Weber-Schaftheitlln  integral 

I  x-X  Ja(px)j^(tx)  dx 

in  which  p  and  t  are  supposed  to  be  positive,  and  where,  to  ensure  convergence 
Re(  a  +  p  +  l)  >  Re(X)  >  -  1  ,  p  /t, 

Re(  a  +  /S  +  l)  >  Re(x)>  0,  p  =  t. 

The  evaluation  of  this  integral  is  a  long  and  complicated  affair,  for  the  details  of 
which  the  reader  is  referred  to  pp. 398-404  of  Watson  (l?44).  We  shall  merely  state 
the  relevant  results. 

It  has  been  shown  by  Sonine  and  Schaftheitlin  that  the  integral 


a+  p  -y 


Ja  -  /?  kx)  Jy  _1(bx)dx 


(3.25) 


has  the  value 


y -a-p&a+PTiy)Ti^fj)  2*1 


(3.26) 


when  b  <  a  and  the  value 


a-p 


2y-„  -P  —*  —  .1r(y  _  a)r(a  ZJTT)  zP'(°’“'y+1i'"/St1C2  ’  (3-27) 

when  b  >  a.  It  will  be  seen  from  these  results  that  the  integral  (3.25)  is  a  function 
of  b/a  which  is  not  analytic  at  the  point  t/a  =  1 .  When  b  =  a,  the  value  of  the 
integral  has  to  be  found  by  a  special  procedure  for  the  details  of  which  the  reader 
is  referred  to  p.402  of  Watson  (1944). 

Some  particular  cases  are  of  special  interest.  If  we  put  a  =  v  +  m  +  1  r 
P  =  1  -  k-m,  y  =  v  +  1,a  =  l,b=p  ,  where  m  is  a  positive  integer,  and  k  and  v 
are  such  that  the  integral  converges  we  have  the  result 


l(i/,m,k,  p  ) 


P  T(y  +  m  *•  l)  _  /  _  _  2. 

r(v  .k)  *(  *’• 1  '  k  -  >’ 


P  >  1, 
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where 


l(  v  ,  m,  k,p  )  =  f  u1  “kJ, 
->» 


P*2m*(u)j>u)da 


Using  the  transformation 


2P1  (  af  P  »  y  ;  z)  =  (l  -  z)  Y  "  “  "^F^  y  -  a  ,y  -  0; y  ;  2) 
(Sneddon,  1956#  p.22)  and  the  definition 

^  (a,  b,  x)  =  F  (-m,  a  +  m;  bj  x) 
m  21 

of  the  Jacobi  polynomial  we  obtain  the  result 


(3.28) 


I(i/,  m,k,  p) 


T  (  v  +  m  +  1 )  0  v  (l 

2k_1  r(w  +  l)r(m  +  k) 


3rm  (k  ♦  V  ,  v  +1 ;  p2),  0  <  p  <  1 


P  >  1 


(3.29) 


To  return  to  Tranter's  series  si  and  sg  defined  by  equations  (3.18),  (3.19)  we 
note  that  it  follows  from  equations  (3.18a)  and  (3.23)  that,  since  m  is  zero  or  a 
positive  integer, 

s  (  v  ,  m,  k,  p  )  =  l(  v  ,  m,  k,  p  )  (3.30) 

where  l(  v  ,  m,  k,  p)  is  defined  by  equation  (3.28)  and  can  be  evaluated  by  means  of  the 
formulae  (3.29).  In  a  similar  way  it  follows  from  equations  (3.19a),  (3.23)  that 

a  (v,«,  n,k)  sja 2  /  x"1  J  (x)j  .  (x)dx 

*  2  J,  ^  +  2m  +  k  y+  2n  +  k  1 

_(H  &_  3ln(k  v )  K  k  v  k  (3.31) 

sr  VfV+tZm+K-t  v  +2n+ 

Another  special  case  of  equation  (3.2 3)  which  is  of  interest  is 

a)  =  /„  Jv(/,x)  Jv-^p(tl)l1”iPdx 

2 

-  -  sin(ipw)K  !  1  (  p  ,  t;  a) 

w  '  v,  p,  p-^p,  — g-  p'  r  *  *  ' 


which,  as  a  result  of  equation  (3.29),  can  be  written  in  terms  of  Heaviside's  unit 
function  in  the  form 
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Sy,  v  ,  v-^p,  -4p^  p  *  t; 


~  gP;(j:r-^P  ^(o-t)  a). 


(3.32) 


In  a  similar  way  we  can  show  that 


* V ,  V  +  1  ,  V •  -  g-p,  — g-  p-1  ^  P  *  a) 


2*4p  tV  -  2P  (p2  _  t2  )2  P 


i  -■  u?  ..  ~*  — .  i-,  -  H(  p  -  t)  -  —  sin(gpjr)K  i  i  (  p ,  t;  a) . 

T(l  +|p)p  +1  V2^  1  w,w  +  1,w^pf-4p-l'^* 


(3.33) 


Further  results  can  be  obtained  by  means  of  the  Hankel  inversion  theorem^  If 
we  apply  this  theorem  to  equations  (3.28)  and  (3.29)  we  can  deduce  that 

/#  p  V+1  (i  -p^aMk+v,  „  +  i  ,p2  )J,(pu)aP  =  -  -rfe-]l-l-1'j-k^u'k  J„ ,2m  +  k 

which,  by  a  trivial  change  of  variable,  we  can  write  in  the  form 


V+X*2  -  xa)kH  3»  (k  +  V  ,  v  +1,  xz/s2)J  (6x)dx 


2  Tiv  +  l)r(m  +  k)  v+  k  _-k  ^  /_  \ 

- TCvVm+\) - L  3  5  Jv  +  2m+  k(Ss>* 


(3.34) 


If  m  =  0  the  Jacobi  polynomial  reduces  to  unity  and  we  obtain  the  simpler  relation 


f*  -  xI)k',J„(«x)dx  =  2kH  r(k)a^t5-kJ„tk(€a). 


(3,35) 


We  also  require  to  obtain  integral  representations  of  two  trigonometrical  series 
involving  Bessel  functions. 

By  integrating  the  function 

cosec(  ifz)eiwz  J  (uz)sin(xz), 
o 

round  the  contour  r  which  consists  of  the  positive  real  axis,  the  positive  imaginary 
axis  and  the  arc,  in  the  first  quadrant,  of  the  circle  |z|  =  R,  with  large  R,  and 

with  indentations  around  the  points  z  =  0,  1 ,  2,  ...,  n,  ...  we  can  easily  show  that 
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e-7ry 

I  (uy)sinh(xy)dy 


oo 

Yt  Jo(nu'l3^nx)  =1  J0 ( tu) sin( tx) dt  -j  3inh(„y) 

n=i  0  0 

Also  it  is  well-known  (Watson,  1944*  p.405)  that 

f 00  _j_ 

/  Jo(tu)sin(xt)dt  =  (x2  -  u2)  2  H(x  -  u) 

j  0 

where  H(x)  is  Heaviside’s  unit  function.  Hence  we  have  relation 

QO 

EJ  (nu)ain(nx)  =  -g(x  --  u)  -  -  /“  ^  M  e’  ^  I  (uy)dy 
>/(x2-u2)  Jo  sinh(wy)  0 

Similarly,  by  integrating  the  function 

cosec(  vz)e1Vz  J  (uz)cos(xz) 

i 

round  the  same  contour  T  we  can  show  that 

PC 

Y  (nu)cos(nx)  =  /  J,  (tu)cos(xt)dt  -  j  e"  vy  I<  (uy)dy 

n  -  .  ®  *  0 


from  which  it  follows  that 


OP 

S'  j  W«»M  »  l  -  -  raa»M  i  (*,)«,. 

nT7  u  'V  (*  -  u  )  Jo  sinh  (iry)  1 


3.3.  Some  Integral  Equations. 

In  solving  dual  series  relations  we  sometimes  have  to  solve  an  integral 
of  the  type 

r  -r  ■h(t)dt  T-a-dx),  a  <  X  <  b 

*  j_f(x)-f(t)J 

vdiere  0  <  a  <.  1  and  f(t)  is  a  strictly  monotonic  increasing  function  in  (a,  b] 
shall  give  a  solution  due  to  Srivastav  (1963c)* 

Consider  the  integral 


X  f'(u)g(x)du  _ 
1  [f(x)  -  f(u)l'-“ 


(3.36) 


(3.37) 


equation 

(3.38) 

•  We 
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If  we  substitute  the  ejqpression  for  g(u)  given  by  equation  (3 *38)  and  interchange  the 
order  of  the  integrations  we  find  that  the  integral  is  equal  to 


h(t)dt 


f  HuJcLll 


t  j^f(u)  -  f(t)J“  jf(x)  -  f(u) 


The  inner  integral  is  easily  shown  to  have  the  value  B(a  ,1  -  a)  =  it  cosech  n  a  „  It 
therefore  follows  that  the  integral  has  the  value 


it  cosec 


/•x 

n  a  h( t)dt 
J  a 


and  hence  that  the  required  solution  is 

h(t)  =  SiSJLiL  ±  [*- 
Ja  1 


f '(u)g(u)du 


it  dt  a  ^f(t)  -  f(u)  1  a 


By  a  similar  method  it  can  be  shown  that  the  integral  equation 


/  r - — ~=g(x)>  a  <  x  <  b 

X  j^f(t)  -  f(x)J 


where  0  <  a  <  1  and  f(t)  is  monotonic  increasing  in  (a,  b)  has  solution 


h(t)  ,  -  sisaa  i  f 

n  dt  J 


* 


f 1 (u Jg(u) 


t  f (u)  -  f(t) 


1  -  a 


(3.39) 


(3.40) 


(3.41 ) 


Two  special  cases  are  of  particular  interest  and  we  shall  consider  these  now. 

If  f(u)  =  -cos  u,  f 1  (u)  =  sinu  and  a  =  ?  and  we  find  that  the  integral  equation 


[  h(t)dt 

a  cos  t  -  cos  x) 


=  g(x)  ,  a  <  X  <  b 


has  solution 


h(t)  -  *  *  J! 


sin  u  gCuJdu 
✓  (cos  u  -  cos  tJ 


,  a  <  t  <  b 


and  the  integral  equation 


h(t)dt 


x  fJ j cos  x  -  cos  t) 


=  g(x),  a  <  x  <  b 


(3.42a) 


(3.42b) 


(3.43a) 


has  solution 


*<*)--£&  /. 


3in  u  g(ujdu 
'(cos  t  -  cos  uj 


a  <  t  <  >>„ 


(3.43*0 


The  classic  case  is  the  one  in  which  f(u)  =  u%  f’(u)  2u.  In  this  instance  we 

find  that  the  integral  equation 


[  x  h(.tj,at _ 

1  a  (xz  -  tz}“ 


:),  0  <  a  <  1  ,  a  <  x  <  b 


has  solution 


h(t)  = 


2  sin 


iin  1 ra  d 
7T  dt  J 


r j  u  gLujfe..  a<  t<  b. 

«  j.  it!  -  «)'-• 


and  that  the  integral  equation 


(3<44&,! 


(3.44b) 


/  ft^v ° g(ll)’  0<“<1’  a<*<b 

J  x  (t  -  x  ; 


(3.45a) 


has  solution 


h(t)  =  -  2 '  a3^-?a  A  f b  0  <  a  <  1 ,  a  <  t  <  h.  (3.45b) 

ir  dt  J  t  (u2  -  t2)1 


Some  results  on  Associated  Legendre  Polynomials. 


In  the  discussion  of  dual  equations  involving  Fourier-Legendre  series  of  the  type 
(2.52),  (2.53)  we  require  certain  properties  of  associated  Legendre  polynomials.  We 
list  these  now. 

Ferrers'  associated  Legendre  function  of  the  first  kind  is  defined  by  the  equation 


Tra  (x)  =  (-0™  (1  -  x2F'"-^-P  (x),  (-1<x<l), 

dx 


(3.46) 


and  n  being  zero  or  positive  integers  and  P  (x)  denoting  the  Legendre  polynomial  of 


degree  r,  and  by  the  relation 


T"”  (x)  -  (-I)""  wr  -r(n-t  II..  Ta 

ra+n  '  '  '  '  T(n  +  2m  +  1 )  m+n  '  ' 


(3.47) 


(UacRobert,  1 947*  PP.125  and  328), 

If  we  assume  that  the  expansion 


f(fi)  =  )  (2n  +  2m  +  l)c  T’ra  (cos  9  ) 

/  1  7  n  m+n  ' 


(3.48) 


is  valid  for  0  s  6  <  ir  and  that  it  can  be  integrated  term  by  term,  then,  using  stan¬ 
dard  integrals  involving  associated  Legendre  polynomials,  we  can  easily  show  that  the 
coefficients  are  given  by  the  formula 


cn  =  ■J(-l)r1'  f  f(x)T^+n  (cos  x)3in  x  dx 
J  0 


(3.49) 


Using  the  results  of  §3.3  we  can  obtain  an  integral  representation  of  the 
associated  Legendre  polynomial.  If  we  write 

_  r  8  R m  (u)du 

T  (cos  e)  =  cot  $6  . . .  0  <  9  <  ir  ,  (3-50) 

J  o  cos  9  -  cos  9  ) 

then  it  follows  from  equations  (3.42a  and  b)  that 

_  1  d  r  u  sin  x  tanm  4xT  (cos  x)dx 

r  ra  (u)  =  — - - — - —  . 

«  du  *o  /ylcos  x  -  cos  uj 

The  integral  on  the  light  is  difficult  to  evaluate  but  it  has  been  shown  by  Collins 
(l96l)  that 


0“)  ■  OsTfM*-  »i)  aln- (Ju)  „o,(iu/— U  aV  safe.* » .>  Th .  (3.5, ) 

nT(n  +2n  +l)  ysin  u  diy  cos^ru 

When  m  =  0,  this  gives  R°n(u)  =  Rn(u)  where 

Rn(u)  J^L.  C03(n  +  i.)u.  (3.52' 

7T 

If  we  insert  this  expression  into  equation  (3.5 0)  we  obtain  the  Mehler-Dirichlet 
integral  for  P^(cos  9),  Similarly,  by  putting  m  =  1  in  equation  (3.50)  we  find  that 


(3.52) 


*'  (u)  =  — 1—  sec(*u)  tan(^u)  ^  3in(n  +  2)u] 

A/2  w  ^  n  +  1  n  ♦  2  3 


(. 3.53 ) 


Rufi-1  Relations  involving  Fourier— Bessel  Series. 

We  shall  begin  by  considering  the  pair  of  dual  series  relations 

00 

Yj  an  Jv^P\)  =  0<P<1, 


(4.1) 
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00 

Ea  J  (pX  )  =  f  (p),  1  <  P  <  a, 

n  y  n  2 


(4.2) 


where  p  and  v  are  real  constants  satisfying  -1  $  p<  1,  v  >  0,  j Xp J  are 
positive  zeros  of  the  Bessel  function  Jy(aX).  The  functions  f  (p)  and  f  (p)  are 
prescribed  and  the  problem  is  to  determine  the  sequence  of  constants  a!  •  This 
pair  of  equations  is  an  immediate  generalization  of  the  pairs  of  equations  (2.7) , 
(2.8);  (2.20),  (2.21 ) . 

We  split  the  solution  of  this  problem  into  two  parts  by  considering  two  special 
problems  the  solutions  to  which  can  be  combined  to  give  the  solution  of  the  general 
problem: 

Problem  (a):  This  is  the  special  problem  which  occurs  when  f  (p)  =  0  in 
which  case  the  equations  reduce  to 


OO 

EX"Pa  J  (pX  )  =  f  (p),  0«  P  <  1 

n  n  n  i  J 

n*  i 

OO 

^an  Jv  (yftXn)  =  0,  1  <  P  «  a. 


(4.3) 


(4.4) 


Problem  (b) :  This  is  the  special  problem  corresponding  to  the  case  in  which 


f  (p)  =0  and  the  equations  reduce  to 

i 


0SP<1 

n*i 

OO 

^  a  n  Jy(pXn)  =  f\,(p),  1  <  p  <  a. 


(4,5) 


(4.6) 


Problem  (a)  was  first  considered  by  Cooke  and  Tranter  C 1 959)  using  a  method 
similar  to  Tranter's  method  of  solution  of  dual  integral  equations.  We  shall  begin 
with  an  account  of  this  method  and  then  give  an  account  of  an  alternative  method  due 
to  Sneddon  and  Srivastav  (l  963)  in  which  the  solution  to  the  problem  is  reduced  to 
that  of  a  Fredholm  integral  equation  of  the  second  kind. 


1 .  The  Cooke -/Tran  ter  Solution  of  Problem  (a 


We  need  not  restrict  v  to  v  *  0  but  we  assume  that  v  is  not  a  negative 
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integer  and  that  v  >  -  1  +  jp. 

If  we  put  k  =  1  -  jp  in  equations  (3.30)  and  (3.29)  we  find  that 


E  - 

n=1 

so  that  if  we  take 


♦  2m  ♦  J  -  jj  (V  m  0>  (,  „  f  ,  a) 


a  =  — - 1 - \ - * -  )  b  J  0  .  t  (\  ) 

n  X1  "  TP  J2  (X  a)  y  +2m+1-§pv  n' 

n  v  +  1 '  n  '  m=o 


(4.7) 


(4.8) 


the  equation  (4.4)  as  automatically  satisfied. 

The  coefficients  have  now  to  be  chosen  such  that  a^  as  given  by  equation 

(4.8J  satisfies  equation  (4o).  Substituting  the  expression  (4.8)  for  a^  into 
equation  (4.3)  and  interchanging  the  order  of  the  summations 

r*  P  0  .  i  (X  )  J  (X  p) 

E  E  1 1)  1 = fi<p)’ (0<  '  ‘  (w) 

m=o  n=1  *n  +1  '^na' 

Now,  by  equation  (3.34),  if  a  is  zero  or  a  positive  integer 


Jy  +  2s  +  1-d-p(Xn)  2zPr(i;  +  s  +  1  )  .1 


r(,  +  i)r(8  +  i-ip);° 


—  f  pw+1(i-p2)4p  X 
i  _\  J  0 


-iP  +  v,v  +  1  ,pz)jy(Xnp)dp,  (4.10) 

where  is  a  Jacobi  polynomial.  Hence  if  we  multiply  equation  (4.9)  by 

'>  V+1  (l  -  p2  )  sP  3a(l  “  +  v  t  v  +  1  »  p  2)»  integrate  with  respect  to  p  from  0  to 

1 ,  and  interchange  the  order  of  integration  and  summation  and  make  use  of  equation 
(4.10),  we  find  that 

00 

bm  Bm(v,  s,  p)  =  B(  w,  s,  p),  (4.11) 

m=o 


where  we  have  written 
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B(  0  ,  3,  p)  =  t1) -  [' 

r(i»  +  i)r(  3+1  -  £p)  Jo 


i/+i 


(i  -  p  2  )”^p  x 


Bm(  ■’.s)p)  = 


-  5P  +  V  ,  v  +  1  ,  p  2)f1(p)dp, 

*^1/  +  2m +  1  -ip  ^  ^ v  +2s  +1  -ip^n^ 


=£ 


n=1 


'2P _ 

X  2  jjj  (X  a) 
n  v  +1  v  n 


(4.12) 


(4.13 


Equation.  (4.11 )  with  s  =0,  1,  2,  3,  ...  provides  a  set  of  algebraic  equations 
for  the  determination  of  the  coefficients  lb  ).  Once  the  values  of  the  b's  have 

f  1  L  m  J  m 

been  found,  the  coefficients  j  j  can  be  evaluated  by  means  of  equation  (4.8). 

The  set  of  equations  (4.1 1 )  can,  under  certain  circumstances,  be  solved  by  an 
iterative  process  as  follows: 

Since  v  is  not  a  negative  integer  and  v  >  -  1  +  jpwe  know  from  equation  (3.31 ) 
that,  for  all  values  of  m  and  s  under  discussion  here,  we  have  the  relation 


(  v  ,  p) 


B  (  v ,  a,  P)  = 


a  m 


where  6 


is  the  Kronecker  delta  end  L 


2p  +  4s+2~p  2v  +  4s+2-p 

(  *  ,  p); 


m.s 


'denotes  the  integral 


(4.14) 


|M 

*  TV+  2s +1  -|p  aj  ‘ 


(4.15) 


Substituting  from  equation  (4.14)  into  equation  (4.11 ),  we  find  that  the  algebraic 
equations  to  determine  the  coefficients  ^  J  become 


b  - 
8 


V  * 

1  mm. 


( v  »  p)  _  £. 


S  2 

a 


(2u  +  4s  +  2  -  p)B(i>  ,  s,  p),  (s  =  0,  1,  2,  ...)  (4.16) 


m=o 


The  iterative  solution  of  this  set  of  linear  equations  (which  can  be  shown  to 
converge  if  a  is  large  enough  -  although  it  is  difficult  to  determine  precise  limits 
for  this)  is  given  by 

00 

\=Yj  tir)  *  (■=0,1,2,...)  (4.1 7) 

r=o 
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where 


b^0'  =  —  (2  v  +  4s  +  2  -  p)E(  u  3  s  ,  p) 


(4-.  18a) 


OO 

b<r>-  r  4">p)  (r  .1,  2,  3,  ...).  (4,13b) 

3  ffi,s  m  - 


The  de t s rmina ti on  of  the  coefficients  s  cf  the  pair  of  dual.  series  relations 

n 

(4.  3)  and  (4,4)  is  contained,  in  equations  (4.8).  (4.,_).  (4.18)  and  (4.1 5) »  in.  2 

practical  problem,  the  chief  difficulty  v/ould  lie  in  the  computation  of  the  integrals 
T"(>,p)  *  ' 


4.2,  The  Ksduct-jon  of  Problem  (a)  ro  an  Integral  Equation. 

We  consider  now  the  method  of  Sneddon  and  Sri^stav  (-1 9 ^3)  by  means  of  which  the 
solution  of  Problem  (a)  can  be  reduced  to  that  of  a  Fredholm  integral  equation  of  the 
second  sc: r,d„  The  cases  u  =  0,  u  >  0  hare  to  be  treated  separately  and,  separate 

consideration  is  required  for  p  >  0  and  p  <  0, 

Problem  (a)  with  0<psl,v  >  0: 

Suppose  tha'  for  0.<  0  <  1  , 


1  J  ( 0  \  )  =  a  ( p) ; 
n  v  n 


(4.19) 


then  the  problem  will  be  solved  if  we  can  determine  cr(p )/ by  the  F ourier-Be s s el 
expansion  theorem  } 


*  a2  J2 


»> 


j  u  a(u)jy(uXn)du, 


(4.20) 


We  now  represent  the  function  a(p )  in  terms  of  an  auxiliary'  function  g(t)  through  the 


r(p)  = 


s(i)(^  -  p"  )  ~  p  d,t 


(4.21) 


and  reduce  the  problem  to  that  of  determining  g(t).  Substituting  from  equation  (4.2l) 
into  equation  (4.20),  interchanging  the  order  of  the  integrations  and  making  use  of 
the  formula  ( 3*3 5)  we  find  that 


21*sPrO  -  5-p)>. ,S,P 


*n  =  a2  J2  (aX  ) 
v+1  v  n  ' 


~  /  tV‘*Pjv  sd)dtt 


(4.22) 
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/< 


Su ss 1  i tuTung  these  values  for  the  coefficients  a  irh'c  squa ton  (4„ 3)  and,  interchanging 
the  order  of  integration  and  simulation  we  get 


/  1  1 

v  -i-D 


r2Fr('  4?)  I  tu-^ug(t)S  „  ±  _2LC  (p,  ";a)at.-=f,(.o),  05,0  <1,  (4.23) 

J  n  9  '  9  *  2iJ3  2  r  ' 


where  S 


5  '=**  9  Ps  "  5  p 


(  p  ^  t;  a)  is  defined 


oy  ecu&f ion  ,  p , '  / )  , 


Maki: 


the  formula  (3-s32)  we  find  that  'his  equation  may  be  written  in  th 


J  r(i  -  r p) 


Vp) 


+t-p 


n(ip^)  r  1 

J  0 


1  <{p»u.a;  du 

h;  us  -  -?P?  *TPV  '  ' 


where  X  is  the  integral  defined  in  §},-!.  Using  the  formula  (3„£)  we  obta.ir 


m  the 


vicn 


s(t)  - 


2- 


r  (1  -  Ip) 


1  U1 

•2>  ~2P  1 


J  -f^rp 


SXU 


jp.  [1"Kv,»,--ip,4p(“,u’a)!  ‘p 


Making  use  of  equations  (3«5)  and  (3»10)  we  see  that  this  equation  reduces  to 


g‘.t)  = 


91 


,P  '  -  V 


- ±1_ 

r  ci  -  Ip)  2 


r] 


+  ~  3  XU  2 
77 


i__  J  +  2P  *  ^  /  ’  v  —in  (  \_ 

sP?  t  /  ^  2‘  g(u)K  1  1 

J0  *  V~*S»  -2p. 


..(t,  u)ftj,  (4.24) 


which  is  a  P.rediiolm  integral  of  the  second  kind  for  the  determination  of  the  functio: 
fi(l')  - 

Problem  (a)  w± th  -1  5  £  <  0,  y  >  0. 

Ijj.  this  case  we  multiply  both  sides  of  equation  t4<,3)  by  p  *  +  and  integrate 
vd.th  respect  to  p  from  0  to  p  to  obtain  the  relation 


L — / 


K  P  '  a  J  (p\  )  =  P  (p)  0  5  p  <  1 

n  n  v  +1  n  'ivr/’ 


(4.23) 
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where 


(p)  =  P~  '  j  n!'+1  (’a)  dr . 

»'  o 


Again  we  assume  that  the  coefficients  J  can  he  represented  by  means  of  the 

equations  (i-.20)  and  (4.2l)  and  hence  by  equation  (4. 22).  Substituting  from 
equation  (4.22)  into  equation  (4.25)  and  interchanging  'be  order  of  summation  and 
integration  we  find  that  equation  (4.2p)  is  equivalent  to  the  relation. 


2"2P  r(l  -  ip)  /  t.y"2Pg(t)8  1  u  iP,  t5  a)dt  =  F  (p) . 

j  n  u  *  '  »  2P»  2?  T 


(4.26) 


Using  equation  (3 .33)  we  find  that  this  reduces  tc 


...i  ir  +1  f«(t)}  p! 

2s2i-J+i(.  J 


a  F.  (p)  ,2  +ip  rl 

— - -  + - sm(  y  p  w)p  /  g(u)  u 

r(i  -  |-p)  w  '  0 


1  „  £p> u;a)du. 
•'  3  V  +  1  J1'*  ZP,  ~2  P  ”  ‘ 


With  the  help  of  equations  (3.5) »  (3,6)  and  (3.10)  -.tc  can  easily  show  that  this 


reduces  ro 


9P  +  1  ,-2i,  -1 

g(t)  =  - - - — - I 

r(i  -  iv)  ° 


>  -b  ~  1  {  /o 


f(u)duj  t 


Jaun(ipJr)  J  g(u)u  2P  t2p  “  "  +1  0(t,  uj  ah 


(4.27) 


which  as  again  a  Fredholm  equation  of  the  second  kind. 

Froblem  (a)  with  -I  <  0,  ^  =  0. 

The  procedure  is  essentially  the  same  as  that  outlined  above  end  the  final 
results  we  obtain  could  be  obtained  merely  by  letting  v  tend  to  zero.  However  a 
slight  modification  of  the  calculations  is  necessary  when  v  =  0  and  the  method  breaks 
down,  f  or  0  <  p  <  1 . 

«  For  -Is  p  <  C  we  assume  as  before  that 


5  J  (p 

r,  o 


(p*n)  =  ~pTp  f  s(t)(te  -  P2  )"2P  dt,  0«P 


and  this  leads  to  the  relation 


•  _  -  - -  /  j0(p  xn)f  r~  [  g(t)(t2  -  p‘)  2  p  at!  ea 

-1  a  J  (aA  )-fc  n43p.>p  J 

1  '  n 
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which,  as  the  result  of  an  integration  by  parts,  can  be  written  in  the  form 


*J2(aX  ) 
i'  n 


["  /  t"pg(t)dt;  -  [  g(  t)dt  f*  J,  tPAn)(ta  ”P  ")  2pdp 

v.  J  o  •  r  Jo 


Using  the  fact  rhat 


I  -P2)  2PdP  =  X°’  t~p  +  pW  2p_1  (-^p)X.2pt  2P  J_lp  (tXn) 

*  O 


we  find,  that 


21  “2px  jp  r(i  -  ip) 

_  _  -  i. _ ~ 

a'J"(aX  'l 


I  J  o  •  •  7 

*  0 


J  ■!  (t  \  )dt 
-i-  '0  n 


(4.28) 


This  expression  for  &n  is  identical  with  the  one  we  should  obtain  by  putting  v  =  0  in 
(4.22).  The  result  of  the  analysis  is  the  same  as  that  given  by  letting  v  — *■  0  in 
equation  (4,27)  so  that  we  obtain  the  integral  equation 


g(t)  = 


2p-  1 


t  r(i  -  Jp  ) 


4,  -iP  - .  [ L  f(u)du!  1 


+  —  3±n(4p  7T )  j  g(u)  u~2  P  i/  +  sPX  -1  1  (t,  Uj  a)du 

n  J0  °>  ~z  P  >  "z  ?»  0 

for  the  determination,  of  the-  function  g(t). 

Problem  (a)  with  p  =  1  ,  v  =  0. 

To  solve  the  dual  series  equations 
00 

)  a.”  1  a  J  (  p  X  )  =  f  (p),  0  <  p  <  1, 


(4.29) 


(4.30) 


OO 

^anJo(^n}  = 


1  <  M  ai 


(4.31) 


we  set 


OO 

Y  a  J  (p  >■  )  =  a(o)  =  -  -  —  [  — 

Lj  "c  ’  pdpjp  4/(t  z-pz) 


(4.32) 
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for  0  $  p  <  1 .  Then  making  use  of  the  formula  for  determining  the  coefficients  in  a 
Fourier-Bessel  expansion  and  using  the  rule  for  integration  by  parts  to  evaluate  the 
relevant  integral  we  find  that 


^  a2J2(a  X  ) 
i  n 


U'i  r i  r*  J  (pX  )dp  ^ 

. h(t)4t  *  x"  I  "M.  *7^} 


The  integral  involving  the  Bessel  function  is  readily  evaluated  and  we  find  that 


^  a2  J2  ( a  X  ) 
i  n 


j  h(t)oos( X n t)dt . 
J  O 


(4.33) 


If  we  now  substitute  this  value  of  a  in  equation  (4.30)  and  interchange  the  order  of 

n 

summation  and  integration  we  obtain  the  relation 


I  t2h(t)s 


_1  _±(  P  ,  t;  a)dt  =  f  (p),  0 


°J  °>  -2»  ~2 


s  p  <  1 


which  because  of  (3.32)  with  v  =  0,  p  =  1  can  be  written  in  the  form 

i. j  ^(t);  pi  =  ~^f  (p)  +  /  th(t)K  i  i  (  P  ,  t;  a)dt. 

~2,2L  J  y?  1  w  Jo  °»  °»  -2»  ~S 

Applying  the  operator  I  i  to  both  sides  of  this  equation  we  obtain  the  Fredholm 

”2»  2 

equation 


.  .  2t  /-tf’(p)dp  4  p 

h(t)  =  j  h(u)K  _j_  _±  ,  (t,  uj  a)dt 

IT  J*^(t  -p)  V2Jo  °»  2,2,1 


(4.34) 


for  the  determination  of  the  function  h(t)  in  terms  of  which  a  can  be  calculated  by 

n 

means  of  equation  (4.33)  and  the  function  a(p)  by  equation  (4.32). 


4.3.  The  Reduction  of  Problem  (b)  to  an  Integral  Equation. 

We  now  turn  our  attention  to  Problem  (b) .  We  begin  by  considering  the  case 
in  which  the  parameter  p  occurring  in  the  dimi  series  equations 


CO 

y  .  j  (,x ) =o, 

L-i  n  n  v  n 


0  $  p  <  1 , 


(4.35) 


^  an  Jv  (  P\ )  =  f00>  1  <  P  <  a. 


(4.36) 
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satisfies  the  inequality  0  <  p  i  1 .  We  shall  assume  that  v  >  0. 
When  0  $  p  <  1 ,  we  make  the  representation 


(  P  \)  =  -  P 


v-1 


-2. 

dp 


f 

JP 


g(t)(tz  -  pz)  sP  dt. 


(4.37) 


Using  the  Pcrurier-Bessel  inversion  theorem  and  performing  an  integration  by  parts 
in  the  integral  involved,  we  find  that 


2 — gPr  (j_:_2P)  X2P  r  g(t)tl'"2Pj  ,  (tx  )at 

n  2t2  \  n  L  y-2P 


a  Jy+1(aXn) 


a2  J2  „  (a  X  ) 
v  +1  n 


[  uf  (u) J  v  (u 
*  1 


\  )du. 
n 


(4.33) 


If  we  substitute  this  expression  for  a  into  equation  (4.35)  and.  interchange  the 
order  of  the  integration  and  summation,  we  obtain  the  relation 

2^pr(l-£p)J  gCth^^S^  v  ,  v  (  p,  t;  a)dt  =  ^1  (p) ,  (4.39) 

where 

$  (p)  =  ~  /  uf(u)S  (p,u;  a)du<: 

From  equations  (4.23  )“  (4.24)  we  see  that  the  equation  (4.39)  is  equivalent  to 
the  Fredholm  integral  equation 

s(t)  ■  -ip  *] 


+  -sin(£p?r)t UV~&  g(u)K  ±  <  (t,  u)du  (4.40) 

V  J  0  v  2  v  2JP»  u 

for  the  determination  of  the  function  g(t). 

The  case  in  which  the  parameter  p  satisfies  the  inequality  -1  ^  p  <  0  can  be 
treated  in  a  similar  fashion.  Multiplying  both  sides  of  equation  (4.35)  by 
p  l*  +  1  and  integrating  with  respect  to  p  from  0  to  p  ,  we  obtain  the  equation 

OO 

Z  X"P"1  ,nJ->ln) 

n-i 


(4.41 ) 
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If  we  substitute  the  expression  (4.38)  for  the  a^s  on  the  left  hand  side  of  equation 
(4,41 )  we  find  that 

24pr(l  -ip)/i1g(t)t,'-^S^v+ny_^)#4p_1(p,  tJ  a)dt  =  *#(p),  (4.42) 

where  0  <  p  <  1  and  ip  ( p)  is  defined  by  the  equation 

*>)  =  -  J  uf(u)Sy  ^  y+1>Vj  _p_1(  P,u;a)du,  0  <  p  <  1. 

Applying  the  analysis  involved  in  the  derivation  of  the  equations  (4.26)-  (4.27)  to 
equation  (4.42)  we  find  immediately  that  the  solution  of  the  problem  is  reduced  to 
determining  the  function  g(t)  from  the  Fredholm  integral  equation 


gP  +  1  +-2  V  -1 

g(t)  =  Z— ■& -  1. 

r(t  -  ip) 


,  -to-i  ['  *] 


^i"(ip.)  I  s(«)p  W-5P  t5 

V  Jo 


-  v+1  K 


a)dU*  (Wf3) 
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5.  Dual  Relations  involving  Dini  Series. 

If  the  function  f(x)  is  defined  in  the  closed  interval  Co,  1]  ,  its  Dini 
expansion  is,  in  general. 


OO 

£ 

m=i 


b  J  (X  x) 
m  m  n 


(5.1) 


where  |^X  1  is  the  sequence  of  positive  roots  (arranged  in  ascending  order  of 
magnitude)  of  the  transcendental  equation 


X  J'(X)  +HJ  (x)  =  0, 


(5.2) 


H  and  v  being  real  constants  with  v  *  -  w.  The  coefficients  b  are  given  by 

m 

the  formula 


2  x: 


"ffl=(x:-  -H2)Jy2(Xm)  -o 


I  tf(t)  J  (tx Jdt. 

j  u  r«i 


(5*3) 


The  expansion  (5.1 )  holds  if  H  +  v  >  0;  if  H  +  p  =  0an  initial  term 

2(v  +l)xvf  ty+1f(t)dt  (5.4) 

Jo 

has  to  be  added  to  the  series,  while  if  H  +  v  <  0  the  equation  (5.2)  has  two 
purely  imaginary  zeros  (+  i  X  ,  say)  and  an  initial  term 


2X2Iy(X  x) 


(X^  +  v2  -H2)I2(Xo) 


1 1  tf(t)ij>(tXo)dt 


(5.5) 


has  to  be  added.  (Watson,  1  % 4;  Chapter  XffTIl) . 
We  now  consider  the  relations 


OO 

Yfn  Xn  =  f>)»  0  $  p  <  a 

n»i 


(5.6) 


OO 

^Cn  Jv(P\)  =f2^)»  «»<  P  <  1. 

n  =  i 


(5.7) 


Solution  in  the  case  -1  £  -p  <  0.  v  >  0. 
We  assume  that  when  0  <  p  <  c 
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oc 

^Cn  Jv(p\  ^  =  “  P  1,-1  ^  S(t)(t2  -pZ)2Pdt. 


n*i 


Prom  equation  (5.3)  it  follows  that 


2\. 


e  = 


"  (\2  -  y  2  +  Hz) Jz(x  )  L  Jc 


y  n' 


j^j*1  tfa(t)  Jv(tXn)dt  +  2^>X^P  r(l  +ip) 


l°g(t)tV+^Jv+h(tXj<Lty 


(5.8) 


If  we  now  substitute  this  value  of  cn  in  equation  (5.6)  and  interchange  the  order 
of  integration  and  summation,  we  obtain  the  relation 

2^pr(l  +*p)  t  v**pg(t)  E'V'  „+ip>ip(p,  t)dt  =  *t(p),  0  <  c  (5.9) 


where  the  function 


^(p)=f(p)+J  tf2(t)Sy,H,  y,  y,  ^dt  <5.10) 


is  known. 

Now  from  equations  (3.24)  and  (3.29) 


♦  f  sm(ip,)  K;>Hi  ^  ..jp.jpto.  t), 

where  H(x)  is  Heaviside’s  unit  function,  so  that  equation  (5.3)  can  be  written  in 
the  form 

which  by  application  of  the  operator  i’1  ±  ,  can  be  written  in  the  farm 

”2.»  ”  2P 
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2-Pt1-u 


v,k>(t'u)dU 


g(t)  - [ii  ,  ff  (p)j  tl  +/-)2P  f  1  uf  (u)S*  H  i  1  (t,  u)du 

r(l  +h>)\j '  j  V  t  /  J  c  2  v+ip,v,&  »  ' 

"  (f)2Pf  8i“(iP')/o  g(u)uV+'PKt,H,  v+ip,v+^,  0  (t’  u)du  * 

Equation  (5.11 )  is  a  Fredholm  equation  of  the  second  kind  by  means  of  which 
we  can  determine  the  function  g(t)  and  hence  determine  the  constants c  J  c 
Solution  in  the  case  0  <  n  $  1  .  v  >  0.  v 

If  0  <  p  «  1  and  v  >  0  the  above  procedure  has  to  be  modified  slightly. 
Multiplying  both  sides  of  equation  (5.6)  by  pv and  integrating  with  respect  to  p 
from  0  to  p ,  we  see  that  equation  (5.6)  is  equivalent  to  the  equation 


oo 

t^f^at,  o$P<c. 

rr.  Jo 


(5.12) 


Substituting  the  expression  (5.8)  for  c  into  this  equation  and  interchanging  the 
order  of  integration  and  summation  we  obtain  the  relation 

2*Pr(l  )j  ty+^PS*jH^+1^+^p^p_1(p,t)g(t)dt=542(p)f  0«p<  c  (5.13) 


where 


(5.14) 


for  0  ^  p  <  c.  Since 


av,  H,  v+1,v+^p,  Ip  -1  ^ 


2iv  tV+^(aZ-  t2)~^P 


r  (i  -  ^p) 


it  follows  that  we  may  write  equation  (5.13)  in  the  form 


’’rffj ^(p) - ; ’ia(ipiy% (5.16) 
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Applying  the  operator  1  ^  i  i  to  both  sides  of  this  equation  and  malting  use 

V  “21 1  1  —  "2 P 

of  the  results  (3.6  ),  (  3.1?  )  we  see  that  this  equation  is  equivalent  to  the 

equation 


^gft)  = 


21",P  &~v 
r(i  +ip) 


,u)du  . 


Substituting  the  expression  (5.14)  for  ^z(p)  we  see  that  this  equation  is 
equivalent  to  the  Fredholm  integral  equation 


01  -  p  . -2  v  -1  f  fP  ,  i 

g(t)=- - 7 -  I0  jl_J  f  uf  (u)du;  t 

r(l  +*p)  °»  2P  Jo  1  J 


?-£p  ti-K-iPfi 


T  -  V  -  g-p  i-  1  , 

— - — — — -  I  uf  (u)  S  „  i  i  (t,  u)du 

+  ip)  Jc  '  ».H,  ^4p,s?pv  * 


r(i  +  ?p) 


21 

Trf  (l  +  -g-p) 


sin (iP#)|o  Kv,  H,y+iE>v*iP,0(t.u)u,^g(u)du  (5.17) 


by  means  of  which  the  function  g(t)  can  be  determined. 
Solution  in  the  case  0<  j  <  1,  v  =  0,  H  >  0. 

We  now  examine  the  dual  series  relations 


^  CnXnP  Jo<  P  0  *  P  <  c, 

n*i 

oo 

^cn  J#(p  An)  =  f2(p),  c<p<1, 

n*i 


(5.18) 


for  0  <  p  ^  1 ,  it  being  assumed  that  the  n  J  3X6  positive  roots  of  the 
equation  (5.2)  with  H  +  v  >  0.  We  make  the  representation 


dt. 


0  <  P  <  c 


(5.19) 
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from  which  it  follows,  "by  equation  (5.3 )»  that 


C"  =  « 


-  ^(1 ')  [/>^(^  -  /0°J>x"fc 

n  0  n  m.  ~ 


=  -[  g(t)tpdt+Xn  f  g(t)dt  f  J  (  pXn)(tz  -  p2)sP  dp 
Jo  Jo  Jo 


f*  j  (A,)(t*  -,*)*» a,  > 

Jo  1  n  X  X1+?P 


•2P  J±_(tx ) 

2  P  n 


so  that 


Cn  (X2  +  h2)jz(  X  ) 

n  o  n 


tf  (t)  J  (tX  )dt 
2'  '  o  n' 


+  22pr(i  +  J°  t*pj^(txn)g(t)dtJ. 


(5.20) 


This  is  exactly  the  expression  we  should  have  obtained  for  cn  if  we  had  put 
v  =  0  in  the  right  hand  side  of  equation  (5.8)  but  it  is  arrived  at  by  a  different 
method;  the  rest  of  the  analysis  follows  the  same  course  as  that  leading  to 
equation  (5.17)  and-  we  find  finally  that  g(t)  is  the  solution  of  the  Fredholm 
integral  equation 


g(t)  =  £ 


1  -  p 


-  -  If,  i  .  j"  [  uf  (u)  du;  tl 

+  i  p)  °»  *P“  1  t  Jo  1  J 


t  r(i  +  £  p) 

24p  ti  -4p  r  i 


i  t  —  2  p  r  •  * 

-  /  uf  (u)  Sn  „  i  -  i  (t,  u) 

+  ip)  Jo  2  °»  H»  ip»°»ipv  * 


r(i  +  ip) 


2  r(iV*P)2P  ‘i"lip  ’Vo  *°>  K-  4p<  4p.  g(u)du- 


(5.21) 
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Solution  in  the  case  v  +  H  =  0. 

The  special  case  in  which  v  +  H  =  0  is  particularly  important.  The  dual 
Dini  series  then  take  the  form 


00 

a°°pV  +  Yj  °n  =  °*P 

n*i 


<  Cl 


(5*22) 


oo 

CQPV  +  ^  C  n  Jy  (pln)  =  f^p),  C  <  P  $  1, 

n*i 


(5.23) 


where  -1  $  p  <  1  and  X  .  X  ,  are  the  positive  zeros  of  J  ,  ,  (X).  The 

12  V  + 1 

constant  a  would  appear  to  he  indeterminate,  hut  we  saw  previously  (Cf.  equations 
(2.27),  (2.28)  above)  that  in  some  physical  problems  its  value  is  known  ah  initio: 
in  other  problems  its  value  may  be  derived,  at  a  later  stage  in  the  analysis , 
from  some  physical  criterion  such  as  that  a  certain  component  of  stress  must 
always  remain  finite. 

When  0  <  p  $  1  the  above  analysis  can  readily  be  adapted  to  meet  the  present 
situation  and  we  find  that  everything  goes  through  easily 5  in  the  end  we  find 
that  the  solution  is  given  by  equations  (5-8)  and  (5-17)  if  v  *4  0  and  by 
equations  (5.20)  and  (5.21 )  if  v  =  0  except  that  in  equations  (5.17)  and  (5.21 ) 
the  functions  K*v  ^  ^  (u,  t)  and  8*  ^  ^  y  ,  g  K  ▼)  are  replaced 

by  Kv  +  1,  /S,  y,  6  (u'  v)  3V  +  1,  p  ,  y,  6  (u’  v)  reapeCtively-  In  other 
words  when  v  /  0  the  equations  have  solution 

c  =  2(  v  +  l)  f  1  tV*  1  f  (t)dt  4  2.r(H.t  £l£liLll2  r  u2t,+pg(u)du  (5.24) 
'c  8  r(v  +  ^p+i)  ■'o 


c  = 


[f1  tf  (t)j  (x  t)at  +  2^x"2Pr(i  +  W  tv+&g(t)j 

KJo  2  v  n  "  Jo  1 

x  (tx^dtj  , 


where  g(t)  is  the  solution  of  the  Fredholm  equation 

g(t)  =  X  (t)  -  f  L(t,  u)g(u)du, 

J  o 

where  the  free  term  x  (t)  is  defined  by  the  equation 


(5.25) 


(5.26) 


b2 


„1  -  p  .-2  v  -  1  r  rP  'j 

*(t)*  r(i.*p>  Xo.iP^[/o  Uf'(u)41’  *J 

P*4p  J  -  y  -  ip  f  1  ,  .  .  . 

+  — - - -  /  uf  (u)S  .  i  i  (t.  u)du 

r(i  +  ip)  Jo  2  v+i,v+ 2P,»»,  2p 
_22_p(y  +  i  )o  r  (^p  +  iv  + 1)  [  i 


( 


l  w  +  i;a  r  i2p  +  iv  +  i /  M 

-  . . . . . —  /  uf  (u)du 

v  +  2) r (^p  +  i) r  (p  +  i  v  )  Jo 


and  the  kernel  is  given  by  the  equation 

21  -^Pt1  -  v  -  ip  uv  +  ip 


L(t,  u)  = 


irT  (l  +  ip) 

22  “  p  r(i»  +  i)r  Up  +  iv  +  i )tp 


sin^pTrjK^  +  1  >  y+  „  + 


-1/-1  2y  +p 

u  a 


r(p  +  i  v )  r(?p  +  v  +  i ) 

In  particular  if  we  put  v  =  0,  p  =  1  ,  f 2(p)  =  0  we  see  that  the  dual 
equations 


a  CQ  +  I  \Cn^o^  V  ~  0  $  p  <  c. 


n*  i 


c  + 
o 


n=i 

where  the  X  n J  are  the  positive  zeros  of  (x)?  have  solution 


=  2  f  ug(u)du, 

J  o 


c  = 


"  X  j2(x 

n  o'  n 


“*  f  g(t)sin(tXn  )dt, 
)  Jo 


where  g(t)  is  the  solution  of  the  Fredholm  equation  (5.26)  with  free  term 


*(t)=£  [* 

v  Jo 


O  ^(t2  -  V  ^ 


and  kernel 


s/z 


2V'  ‘  r°°  K  (y) 

L(t,  u)  =  2  «u  +1  —  I  /  "  sinh  (ty)sinh  (uy)dy. 

\  w/  Jo  I^y) 


(5.27) 


(5.28) 

series 

(5.29) 


(5.30) 


(5.31) 


(5.32) 


(5.33) 
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Solution  in  the  case  p  =  -d  .  v  =  H  =  0. 


It  is  difficult  to  derive  the.  solution  in  the  general  case  when  H  +  y  -  0 
and  -1  s  p  <  0j  we  shall  consider  only  the  two  cases  which  are  of  most  interest 
physically.  The  first  problem  we  shall  consider  is  the  solution  of  the  dual 
equations 


ac  X  "1c  J  (pX  )  =  f(p),  0  $ 

o  / .  n  non  r" 


p  <  c 


(5.34) 


^  °n  =  °» 


B  <p  (1 


(5.35) 


in  which  the  function  f(p)  is  prescribed  in  (0f  c)  >  a  is  a  given  constant  and 
X  i  ,  Xa,  X3,  ...  are  the  positive 


J  (X)  =  0. 


We  suppose  that  when  0  «  p  <  c 


c  +  )  c  j(pO=-U 


no  n' 


tcltidt 


p  3p  JR  y(tz  -  p2) 


(5.36) 


=2  \  g(t)dt, 

*  O 


C"  =  JZ(X  ) 
o'  n ' 


/  c  g(t) 
J  o 


cos(tX  )dt. 
n 


(5.37) 


Substituting  these  expressions  for  the  coefficients  cn  into  the  left-hand  side  of 
equation  (5.34)  we  find  that 


C  C  /  oo 

2a  J  g(t)dt  +  J  g(t)dtf  2  ^ 

V  n«i 


J  (pX  )cos(tX  )'\ 

— - 7*-  )  =  f(p)>  0  <  p  <  o  .  (5.38) 


However  it  can  be  shown  that 

„  r  j#(p\)c °B(tx  ) 


i  i 

- - - -  =  (p*  -  t2)  2  H(p  -  t)  -  2(1  -  t*)* 

V  X  / 

a  ’  n  ' 


X  X 
n  o'  n 


2  /*  K,(y)  f  -) 

o  00811  *y [2I,Cy)  -  yi0(py)j  ay  (5.39) 
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oo 

P  *  Tj  °"  Ji(pXn)  =  °»  C<P  «  1 


(5.46) 


where,  as  before,  the  constant  a  and  the  function  f(p)  are  prescribed,  but  now 
X  1  ,  X  2  ,  ...  are  the  positive  roots  of  the  equation 


J  (X)  =  0. 
2 


We  make  the  assumption  that 


coP  + 


f  cj  (px  )  =  -  — ■  r  -&^dt2  ,  0* 

Z-/  "  1  "  3  p  ip  V(t  -P  ) 


p  <  c 


(5.47) 


from  which  it  follows  that 


=  8  f  ug(u)du, 
J  o 


°n  =  J2(X  )  Jo 

i  n 


fC 

/  g(u)sin(x  u)du, 

Jo  n 


(5.48) 


Substituting  these  values  for  the  coefficients  a^  in  the  left-hand  side  of 


equation  (5«33)  we  obtain  the  relation 


c  c  °° 

8  ap  j  ug(u)dn  +  J  g(t)dt  fz 

n«i 


]  =  fw,  o « „  <  (5.w) 


Considering  the  Dini  expansion  of  the  function  (of  p  )  defined  by  the  right  hand 
aide  of  equation  (5.50)  below  it  may  be  established  that 

ZJ  (pX  )sin(t  X  )  tH(p  -  t) 

-4 - - 0-  = - 4tpV(l  -  t2) 

,  x„J2(x„)  P/vV-t2) 

2  /  »  K  (y)  .  r  . 


2  /  ”  h  (,y;  r  -n 

"  “  /  ~~  sinh(ty)  4pl  (y)  -  yl  (py)  d;y.  (5.50) 

y  4  2  1  J 


n  J  o  yl^  (y)  L  2  1  J 

Inserting  this  expression  for  the  infinite  series  occurring  in  the  left-hand 
side  of  equation  (5.37)  we  see  that  this  latter  equation  reduces  to 


/  ■  =  pf(p)  +  f  g(u)  H(u,  p  )du 

*  o  /(p2  -  t2)  Jo 


where 


H(u,.p)  =  4up2V(l  -  u2)  -  8ap 2  u  +  -  -^iinh(uy)^4p2I2(y)-pjii:i.(fly)J  dy. 
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This  in  turn  leads  to  the  Fredholm  integral  equation 


where 


r  c 

tg(t)  -  /  K(u,  t)du  =  x  (t) 
J  o 

Ht) 

7T  dt  *  0  /l/(  t  -i 


V(t  -pz) 


(5.51) 


(5.52) 


(u,  t)  =-± r 

ir  dt  J  0  /y/(tz  -  p2  ) 


=  ^«t!V(t  -  u!) 

7T 

4  r°°  K2(y)  r  2  -n 

+  “  /  - 7T  8t  I.(y)  “  ty  3inh  (ty)  sinh  (uy)dy. 

»  Jp  yi  (y)  4  J 


Using  the  fact  that 


I  -  <)  ay  -  .  -  <, VO  -  •> ) 


we  find,  that 


.v  ,w  z,  4  /•«  K  (y)  p 

u,  t;  =  ut  Cl  -  2  a)  +  —•  - 3ut  I  (y)  -  ty  sinh(ty)sinh(uy)  dy  (5.53) 

*  *  Jo  yl„(y)  L  J 


0  yiz(y) 


6.  Dual  Relations  involving  Sin  Series. 

We  saw  in  %2  (Cf.  equations  (2.50)  and  2.51)  that  certain  mixed  boundary 
value  problems  concerning  plane  harmonic  functions  lead  to  the  solution  of  dual 
series  relations  of  the  type 

oO 

^  '  n*3  an  sin  nx  =  f(x),  0  (  x  <c  (6.1) 


a  sin  nx  =  0, 
n 


c  <  X  ■£  TT. 
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In  this  section  to  shall  consider  only  the  cases  p  =  +  1  since  these  are  the 
ones  which  occur  most  frequently  in  applications. 

6.1 .  The  case  in  which  o  =  -1  . 

If  to  assume  that  when  0  ^  x  <  c 

(6.3) 


O0 

l 

n  =  i 


sin  nx 


=  g(x) 


then  the  constants  a  whose  values  we  wish  to  determine  are  the  Fourier  sine 
n 

coefficients  of  the  function  defined,  by  the  right-hand,  sides  of  equations  (6.2) 
and  (6.3).  Hence  we  have 


.  .4  r 

n  _  / 

V  Jq 


g(t)  sin(nt)dt. 


(6.4) 


The  problem  of  determining  the  an  (or  the  g(t)  Beems  to  have  been  considered 
first  by  Tranter  (l 959  b)  who  made  the  integral  representation 


g(t)  =  e(i)  sin(^-t)  (i  -  sin2  cosec2 


f  1  £'(s)  3in(?  x)ds 

JT  V(s2  -  Tz) 


(6.5) 


where  T,  the  loTOr  limit  in  the  integral  on  the  right  is  an  abbreviation  for 
sin^- 1  coseo^-c.  Substituting  this  expression  for  g(t)  into  equation  (6,4)  and 
then  substituting  the  resulting  expression  for  an  into  equation  (6.1 )  it  is 
possible  to  find  an  integral  equation  for  £  (s)  whose  solution  enables  us  to 
determine  the  function  g(t)  and  the  set  of  constants  ^ a  n  J  .  Tranter's  analysis 

is  very  complicated  and  for  that  reason  we  shall  not  repeat  it  here.  It  turns 
out  in  the  end  that 


5(a) 


.  2.24 

1  -  s  sin 


it s.sin2 


(6.6) 


The  function  £  (s)  is  determined  in  terms  of  the  known  function  f(x)  by  means  of 
this  equation;  g(t)  is  then  found  from  equation  (6.5)  and  the  constants  a  from 

(6.4). 

% 

Recently  Williams  (l  963)  has  given  a  much  more  direct  and  simple  method  of 
solving  the  pair  of  equations  (6.1 )  and  (6.2).  Substituting  from  equation  (6.4) 
into  equation  (6.1 )  we  find  that  the  unknown  function  g(t)  is  the  solution  of  the 
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integral  equation 


f  c 

/  g(t)  K(t,  x)dt  =  irf(x),  0  <  x  <  c 
J  o 

in  which  the  kernel  K(t,  x)  is  defined  by  the  equation 

OO 

K(t,  x)=2£  n  1  sin(nx)  sin(nt) 


(3.7) 


=  log 


tin  +  t) 
lin  -g-(x  -  t) 


(3.8) 


By  an  ingenious  method  Williams  reduces  this  kernel  to  a  form  which  at  first 
sight  seems  to  make  the  integral  equation  (6.7)  more  complicated  Jjut  which,  in 
fact,  enables  us  to  solve  it  in  explicit  form.  We  consider  the  integral 

.  min(t,  x)  tan(-^  u)du 

X  =  /  -  . 

JO  _  _ j_\  / _  _ 


>\/^(coS  u  -  cos  t)(cos  u  -  cos  x)J 


The  result  will  be  symmetrical  in  x  and  t  so  for  convenience  we  assume  that 


x  >  t  and  consider 


tan  £  u  du 


>^j(oos  u  -  cos  t)(oos  u  -  cos  x) 


If,  in  this  integral,  we  change  the  variable  of  integration  from  u  to  v  where 


y2  =  tan2  £  t  -  "tan2  £  u 


find  that 


I  =  sec(£x)  sec(£t)  f 

J  o 


tan-£  t 


^(v2  ♦  tan2  £x  -  tan2  £t  ) 


This  integration  is  elementary  and  leads  to  the  result 


I  =  £  sec(-g-x)  sec  (£t)  log 


tan  ■Jrx  +  tan  £t 
tan  -  tan  £t 


=  £  sec(£x)  seo(£t)  K(t,  x). 


In  other  words 
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f  min(  t ,  x)  tail  2”U  an 

K(t,  x)  =  2  coa-g  x  cosi t  /  ■  11  ”  .  (6.9) 

*  0  (°osu  -  aost)(cosu  -  cos  x)  j 

If  we  now  substitute  this  expression  into  equation  (6.7)  and  interchange  the  order 
of  the  integrations  in  the  integral  on  the  right  hand  side  we  find  that  equation 
(6.7)  can  be  written  in  the  form 

/•*  tan(-g-u)du  r  c  cos(^t)  g(t)dt 


tan  |u  du 


.  ix[ 

*  Q 


Q  /\/(cos  u  -  cos  x)  J  u  V(  cos  u  -  cos  t) 


/•  c  cos  (it)  g(t)dt  ^  , 

/  -  — . — - —  ■■  ■  -  -iirf(x),  0  $  x  <  c.  (6.i o) 

J u  a/(cob  u  -  cos  t) 


(Cf.  Pig.  5  below  ).  Hence  if  we  write 


(u)  =  r 

J  u 


c  cos  (£t)  g(t)dt 
1  d( cos  u  -  cos  t) 


equation  (6.1 0)  becomes 


tan(-g-u)  (S-(u)du 

/  —  ■  . .  1  =  i  IT  sec(ix)  f(x),  0  $  x  <  c 

i0  V (<JOS  u  -  cos  x) 


(6.11 ) 


(6.12) 


Fig.  5 

Using  equations  (3.42a  and  b)  we  see  that  the  solution  of  equation  (6.12)  is 

d  f  u  sin(ix)  f(x)dx 

G(u)  =  cot  i  u  —  j  -  (6.1?) 

du  i  o  V (cos  x  -  cos  u) 

and  using  equations  (3. a  and  b)  we  see  that  the  solution  of  equation  (6.11 )  is 


_  -  sec  it  d  r°  sin  u  G(u)du 

ir  dt  J  t  \^(cos  t  -  cos  u) 


(6.14) 


We  can  write  this  result  in  another  way.  Integrating  with  respect  to  t  from  x 


50 


to  c  we  find  that 


r  c  sec  ix  rc  sin  u  Cr(u)du 

/  «<*>*- - /  -T. - 7 

i  x  ir  *  x  vCcos  x  -  cos  u; 


Arc 

—  +  —  /  sinuGr(u)du  x 
u)  27rJ  x 


r  u  sec  gy  tan  gy  dy 

x  /  -  . 

J*  ✓(  cos  y  -  cos  u) 

The  y-integral  is  easily  reduced  to  the  form 

1  ,u  sec2  iy  tan  £y  dy  _r  2.  ,  .  ai  2l  T|y=u 

V^2  J  x  ^(l  -  cos2  gu  sec2  gy)  L 


=  sec(gx)  sec2  (gu)  >\/( 


cos  x  -  cos 


(6.15) 


Inserting  this  expression  into  the  double  integral  on  the  right-hand  side  of 
equation  (6.15)  we  find  that 

rc  2  rc  tan(gu)  &(u) 


rc  2  rc  tanlgujGluJ 

/  g(t)dt  =  -  cos(gx)  /  — - - - 

*  x  ir  *  x  A/(cos  x  -  cos  u 


Differentiating  both  sides  of  this  equation  with  respect  to  x  we  find  that 


2  d  r 

g(x)  = - cos(g  x)  1 

ir  dx  J : 


c  tan(gu)  G-(u) 


l  cos  X  -  cos 


-du,  0  £  x  <  c. 


(6, 46) 


where  G-(u)  is  defined  by  equation  (6.13). 

An  alternative  solution  has  been  derived  by  Srivastav  (1963)  based  on  the 
integral  representation 


d  t  [' 

=  -  — ■  COS  g  X  / 

dx  J  x 


,  0  <  x  <  c 


(6.17) 


L  COS  X  -  COS 


of  the  function  g(x)  introduced  in  equation  (6,3).  Substituting  this  form  into 
the  right-hand  side  of  equation  (6.4)  and  integrating  by  parts  we  find  that 


in  =—  j  cos(nx)  cos(gx)dx  f 

ir  J  o  J : 


x  ^(cos  X  -  cos  t) 


Interchanging  the  order  of  the  integrations  we  find  that 

2n  rc  r  t  cos(g  x)  cos(nx) 

an  =  —  /  (t)dt  /  -  dx  . 

w  « o  Jo  ^(cos  x  -  cos  t) 


(6.18) 
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Prom  the  representations 


.  .  ^ r  cos(n  -  ?)x  dx  4/2  fv  sin(n  -  -§-)x  dx 

P  (cos  t)  =—  /  ■  ■  -  I  ——————  (6 

n'"1  it  *  #  4/(003  x  -  cos  t)  v  ^/(cos  t  -  cos  x) 

of  the  Legendre  polynomial  (Whittaker  and  Watson,  1927,  p.315)  we  immediately 


(6,19) 


deduce  that 


4/2  /•  t  cosCg-x)  cos  (n  x)  t/2.  rir  cos(^-x)  sin(  n  x) 

—  / -  dx  = -  /  - 

v  Jo  V(oos  X  -  COS  t)  IT  it  V^(gOS  X  -  cos  t) 


=  i(p  (t)+P(t)"]  (6.20) 

L"-1  n  J 


so  that 


‘.=>^  j"  *<*>[*„-,  <*>  dt- 


(6.21) 


Substituting  this  expression  for  into  equation  (6.1 )  we  find  that  this 
equation  is  equivalent  to 

f  <^(t)  L(t,  x)dt  =  f(x)*  0  <  x  <  c 


(6.22) 


where  the  kernel  L(  t ,  x)  is  defined  by  the  equation 

OO 

L(I’ ±  I  ^Pn(cos  t)  ♦^n_1  (cos  t)^  sin(nx). 
n»i 

Now  it  follows  immediately  from  the  second  of  the  equations  (6.20)  that  the 

series  on  the  right  hand  side  of  this  equation  is  the  Fourier  half -range  sine 

series  of  the  function  /-1  \  TT/-  .  \ 

cosVg-x)  H^x  -  t) 

L(t,  x)  =  -  0  $  x  $  v. 

\/(cos  t  -  cos  x) 

Substituting  this  form  for  the  kernel  L(t,  x)  into  equation  (6.22)  we  see  that  it 
is  equivalent  to  the  equation 


fx  *.Ct)at 

J 0  /SCcos  t  - 


cos  x) 


=  sec(^x)  f(x),  0  s  x  <  c. 


Using  equations  (3.42 a  and  b)  we  can  invert  this  equation  to  obtain  the  formula 


2  d  r  t  sin(^-x)  f(x) 

<£(t)  = . .  . .  ■  ■  1  dx 

w  dt  J  0  a/(oob  x  -  cos  t) 


w  dt  J  0  /y/(cos  x  -  cos  t) 

for  the  determination  of  the  function  <f>  ( t).  If  we  introduce  the  function  G- 
defined  by  equation  (6.1 3)  we  find  that 

$(t)  =  tan  £  t  Cr( t) . 


(6.23) 


(6.24) 
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Substituting  this  expression  for  <f>  (t)  into  the  right  hand  side  of  equation  (6.17) 
we  obtain,  once  again,  the  equation  (6.1 6)  for  the  function  g(x)  in  the  interval 
(o,  c). 

6.2.  The  oase  in  which  -p  =  +1 . 


To  solve  the  equations 


CO 

y'  n  &n  sin  nx  =  f(x),  0  <  x  < 


CD 


sin  nx  =  0, 


C  <  X  <  V 


(6.25) 


(6.26) 


we  integrate  both  sides  of  equation  (6.25)  with  respect  to  x  to  obtain  the 
equivalent  relation 


OO 

^  a  ^  (l  -  cos  nx)  =  j  f(u)du,  0 


<  X  <  c. 


(6.27) 


If  we  now  assume  that,  in  the  interval  0  <  x  <  c, 


oo 

y  a  sin  nx  =  sin(-g-x)  f  g(t)dt 

L—1  °  J  x  a/(cos  x  - 


cos  t) 


(6.28) 


it  follows  from  the  theory  of  Fourier  series  that 


j-n  =  —  J  sin(nx)  sin(-g-x)  dx  J 

Ire  ft  cos(n  - 

=  -/  g(t)dt  /  - j 

nJ  O  Jo  A/C 


J  x  /\/(cos  X  -  cos  t) 
t  cos(n  -  5-)x  -  cos(n  +  ^)x 


(cos  x  -  cos  t) 


Making  use  of  the  integral  representation  (6.19)  for  the  Legendre  polynomial  we 
see  that  this  equation  can  be  written  in  the  form 


"  //  2 


j  g(t)  jy,.,  (cos  ^  “  pn  (cos  tj^dt. 


(6.29) 


If  we  substitute  this  expression  for  the  coefficients  into  equation  (6.27)  and 
then  interchange  the  order  of  summation  and  integration  we  obtain  the  integral 
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equation 


rc  rx 

J  g(t)  S(x,  t)dt  =  J  f(u)du,  0 


<  X  <  c 


where 


°°  -s 

s(x’ 4)  *  ±  I (1  -  cos  nx)  j^P^  (cos  t)  -  Pn  (cos  t)J 
n  « i 


00 

=  JT+Yj  [Pn  (C0S  ^  "  Pn  “  1  (OOS  ^  ] 


cos  nx,. 


(6.30) 


By  considering  the  half -range  Fourier  cosine  series  of  the  function 

sin(£  x)  H(x  t-  t) 

-  (6.31) 

^/(cos  t  -  cos  x) 

we  can  show  that  S(x,  t)  is  equal  to  this  function.  Substituting  the  expression 
(6.31 )  for  S(x,  t)  into  the  left-hand  side  of  equation  (6.30)  we  obtain  the 
integral  equation 


t)  dt 


(cos  t  -  cos  x) 


=  cosec 


(•?  x)  f  f (u)du 

J  O 


(6.3 2) 


for  the  determination  of  the  function  g(t).  The  solution  of  this  equation  is 
given  by  equations  (3.42a  and  b)  in  the  form 

2d  ft  cos(|-  x)  r  i 

g(t)  = - I  -  -  - - — .  dx/  f(u)du.  (6. 

V  dt  )  o  V( cos  x  -  cos  t)  Jo 


(6.33) 


6.5.  Solution  of  Equations  of  the  Second  Type. 
The  solution  of  the  equations 


00 

P  an  sin  nx  =0,  0 


<  x  <  o 


(6.34) 


op 

^  an  sin  nx  =  f(x),  o  <  x  *  1 


(6.35) 


can  readily  be  deduced  from  that  of  equations  (6.1 )  and  (6.2).  If  we  change  the 
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(\  n  ♦  i  p 

“1  J  n  an  =  ‘n* 

y  =  v  -  c  we  find  that  these  equations  are  equivalent  to 


T'  n”p  A  sin  ny  =  f  (y), 

U  "  1 

n«i 


r 


L  n 

n»i 


A  sin  ny  =  0, 


0  $  y  <  y 


y  <  y  <  ir 


(6.36) 


(6.37 


where 


f  (y)  =  f(  ir  -  y).  (6.38) 

If  p  =  1  the  solution  of  equations  (6.36)  and  (6. 37)  is  given  by  equations 
(6.2l)  and  (6.23)  in  the  form 

A  =  —  [  <f>(t)  f  P  (cos  t)  +  Pn  (cos  t)l  dt 

n  4/2  J  a  L 


(6.39) 


where 


2  d  r 

*(t)  =  - 

*  dt  J  1 


2  d  ft  sin(£  y)  f  (y) 


v  dt'*  ^(cos  y  -  cos  t) 


dy. 


(6.40) 


Reverting  to  the  original  variables  we  find  that 

rjr  r  ^ 

j  0(t)  |^Pn  (cos  t)  -  Pn  _  ^cos  t)jdt 


with 


n  V  2 

2  d  f(x)  cos(|-x) 

^(t)  = - - - -  dx 

jr  dt  ^  t  v  (°os  *  -  cos 

is  the  solution  of  the  pair  of  dual  integral  equations 


(6.4l) 


(6.42) 


<0 

n  an  sin  nx=  0, 
n*  i 

oo 

^  a  n  sin  nx  =  f(x), 
n«i 


0  $  x  <  c 


o  <  x  «  v. 


(6.43) 


(6,44) 


If  p  =  -4  the  solution  of  equations  (6.36)  and  (6.37)  is  given  by  equations 
(6.29)  and  (6.33)  in  the  form  ^ 

An  =^2  /  *(t)[Pn-,(cos  t)  “  Pn(oos  t)]dt 
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where 


z  a  f 

=  -  — 

v  dt  J  i 


2  d  ft  cos(-g-y) 


dy  f  f  (u)du 
J  O 


7T  dt  J  o  ^(cob  y  -  cos  t) 

Reverting  to  the  original  variables  we  find  that 

a  =  — [  ift( t)  [  P  (cos  t)  +  P  (cos  t)j  dt 

n  /1/2  ic  l  "  n-  1  J 


(6.45) 


with 


2  &  t  sin(^-  x)dx  .7r  sin(^- t).ircos(^u)f(u)du 

^(t)  = - —  /  . . .  /  f(u)du  = - '  -j  — -  (6.46) 

jt  dt  J  t  /V^(  cos  t  -  cos  x)  Jx  ir  i 


rf(cos  t  -  cos  x) 
is  the  solution  of  the  pair  of  dual  integral  equations 
sin  nx  =  0,  Q  $  x  <  c, 


.  ^(cost-cosu  ) 


)  n  a  su 

L 

n«  i 

00 

a  sin  nx  =  f(x) 


C  <  X  «  V. 


(6.47) 


(6.W) 


n*i 


6.4.  Sine  Series  Analogous  to  Dini  Series. 

Finally  we  shall  consider  the  solution  of  the  pair  of  fhml  integral 
equations 


OO 

^  (n  -  ?)P  a  n  sin(n  -  *)  x  =  f  (x)  ,  0 

n*i 

oo» 

Yj*-  n  sin(n  -  £)  x  =  tg  (x)  , 


s  x  <  c, 


C  <  X  $  V  . 


(6.49) 


(6.50) 


n«i 


where  p  =  +  1,  following  the  method  of  Srivastav  (1965b). 

.Case  (a):  p  =  1  ,  fo  (x)  =  0.  To  find  the  solution  in  this  case  we  n««nn»»  that 


when  0  <  x  <  c 


n»i 

from  which  it  follows  that 
2  f  c 


Zr  c  g  ( t)  d  t 

an  sin(n  -  £)x  =  sin  x  /  — * - -,0<x<c  (6.5l) 

*  x  V(cos  x  -  cos  t) 


=  -  J  ain(n  -  £)x  .  sin  x  d x  j 


c  g^tjdt 
x  /^(cos  x  -  00s  t) 
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--  t)dt  [ 

IT  J  O  J  C 


t  coa(n  -  — )x  -  cos(n  +  ^-)x 
2 _ _ 

j  h/(ooa  x  -  cos  t) 


Using  equation  (6. 1 9)  we  find  that 


^=-$2  j  «,(*)  |j  “  p  (coa  *)]dt 

^  1°  g^t)  |pn_2  (cos  t)  -  Pn  (cos  ‘1  dt,  n  *  2  . 


(6.52) 


If  we  integrate  both  sides  of  equation  (6.4-9)  from  0  to  x  <  c  we  get  the  relation 


^  a^p  -  cos(n  -  £)xj  =  j  f i  (u)du. 


(6.53) 


Substituting  the  expressions  (6.52)  for  the  coefficients  an  and  interchanging  the 
order  of  summation  and  integration  we  obtain  the  relation 


&  /;  g,(t)  s(t,  x)dt  =  /  *  f  1(u)du 


(6.54) 


where 


w 

S(t,  x)  =  p.  -  cos(£x^  p  -  Pi(coat)J+^  |^Pn_2(cost)-Pn(cos  t)jp-cos(n-4)xJ  . 

w  n  *  2 


Now  it  is  easily  shown  that 


V2  sin  x  Hlx  - 
/ V^(cos  t  -  cos 


00 

=p  +  Pi(coa^|cos(2-x)+^pn(cost) -Pn.|cos  t)Jcos(n-^)3 

x)  1-  n>2 


and  that 


OO 

E  1%  ^(cos  t)  -  P^COS  t)J  =  1  +  P^cos  t) 


from  which  it  follows  that 


S(t,  x)  =  2^1  -  cos(^xfj 


W2. sin  x  Hix  -  t^ 

+  Vicos  t  -  cos  x)  * 


If  we  substitute  this  expression  for  S(t,  x)  into  equation  (6.54)  we  obtain  the 
integral  equation 
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,lnx  l 


x  g  (t)dt 


o  V(  cos  t  -  cos  x) 


=  J  (u)du  -  4/2  J°(l  -  cos-g  x)g^  (t)dt 


for  the  determination  of  g  (t).  If  we  write 

1 


yi  =  [  g,(u)du 

J  o 


(6.55) 


and  solve  this  integral  equation  by  making  use  of  equation  (3.42b)  we  find  that 


1  a  s 

5i(t)  = - 

7T  dt.  J 


0  VXcos  x  *“  cost 


/> 


d  r  t  1  -  COS  -g-X 


)du - 1 - , 


u  dt/  0  f^(  cos  x  -  cos  t) 


■d  x.  (6.5 6) 


To  obtain  the  value  of  the  constant  we  integrate  both  sides  of  this 
equation  with  respect  to  t  from  0  to  c,  In  this  my  we  get  the  equation 


r  V2 

r  c  (l  -  cos  g-  x)dx  -> 

1  1 

I  1 

f°  ix 

rx 

i  ,  J 

O  V(c°S  X  -  cos  c)  J 

rJ 

0  ^(cosx  -  cosc)  J 

0 

(6.57) 


from  which  to  calculate  y  . 

Case  (b):  p  =  1,  f^(x)  =  0. 

In  a  similar  way  we  obtain  the  solution  in  this  case  by  writing 
f-1 .  d  f  x  g  (t)  dt 


L(n  -  -g)a  sin(n  -  -g)x  =  -  —  sin  x  / 

"  dx  J, 


«  V(  cos  t  -  cos  x) 


,  c  <  x  <  V,  (6.58) 


which  is  equivalent  to  assuming 


r? 

/ 

J  0 

|  1  +  P  1  (cos  t)J  dt  "|^ 

/•w 

/.  8'(t) 

j*n  (cos  t)  -  Pn_2(cos  t)J  dt. 

(6.59) 


By  a  procedure  similar  to  that  outlined  for  case  (a)  we  can  show  that 

1  d  fV  f  (x)  dx 

g2(t)=-~  /  - * - . 

ir  dt  J  t  /y^(cos  t  -  cos  x) 

Case  (c):  p  =  -1 r  f  (x)  =  0. 


i  a  f 

,(*)  - - 

1  it  dt  Ji 


(6.60) 


To  derive  the  solution  in  this  case  we  begin  with  the  assumption  that  when 
0  $  x  <  o 

Zd  r  c  g  (t)  d  t 

a  sin(n  -  g)x  - - /  ■'  . .  -  (6.6l ) 

n  d xJx  V(cos  x-  cos  t) 


1 
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i, 


7 


I 


which  is  the  same  am  assuming  that 


CL  = 

n 


n  -  p  c 

g.(t)  p  (  cos  t)dt 

Jo  1 


A/2 


n-, 


(6.62) 


(n  =  1  ,  2_,  . ..).  Substituting  this  expression  for  into  the  series  on  the 
left-hand  side  of  equation  (6.49)  with  p  =  -1  and  using  a  procedure  similar  to 
that  adopted  for  case  (a)  we  can  show  that 

1  d  r  t  sin  x  f  (x)  d  x 


1  a  f 

:,(*)=--/ 

tt  dt  Jo 


.  V(  cos  t  -  cos  x) 


(6.63) 


Case  (d):  p  =  -1,f(x)=0» 


To  solve  the  equations 

(n  -  ±)  '  a n  sin(n  -  £)x  =  0,  0  5  x  <  c. 


l 

n*  t 


00 

^  a  nsin(n  -  £)x  =  f2(x),  c  <  x  s  v . 


(6.64) 


(6.65) 


we  integrate  equation  (6.65)  with  respect  to  x  from  x(  >  c)  to  ir  to  obtain  the 
equivalent  relation 


L (ri  ~ 


^rx  -  g;  ap  cos 


(n  -  £)x  =  f  f  g(u)du. 

*  TT 


(6^ 66) 


We  now  set 


<P 


X  g2(t)  dt 


c  /v/(cos  t  -  cos  x) 


c  <  X  <  ff 


(6.67) 


which  is  equivalent  to  the  assumption 

f  g2(t)  PR_1  (cos  t)  d  t. 
j  c 


n  -  i 


V2 


(6.68) 


Inserting  this  expression  for  the  coefficient  a  R in  the  series  on  the  left  hand 
side  of  equation  (6.66),  interchansing  the  order  of  integration  and  summation 
and  making  U3e  of  the  result 
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H(t  -  x) 

\/(cos  X  -  cos  t) 


ZP  (cost)  cos(n-^)x 
n  “  i 


we  find  that 

g2(t) 


1  A. 

7T  d  t 


sin,  x  dx 
tf(cos  t  -  cos  x) 


After  a  little  manipulation  it  can  be  shown  that 

m  f  (u)  du 

_ a v  _ _ 

t  ^/(cos  t  -  cos  u) 


g  (t)  =  - 

2 


2  sin  t 


(5.69) 


(6.70) 


(6.71) 


7.  Dual  Relations  involving  Cosine  Series. 

In  this  section  we  shall  consider  dual  series  relations  of  the  type 

oo 

c*-  ^ 

^  X.  ao  +  ^  a„  cos  nx  =  f  (x),  0  $  x  <  c 

n  =  t 

OO 

^  a  ncos  nx=  fa(x). 


2  a0  + 


C  <  X  <  K 


n~'> 


and  other  dual  equations  involving  cosine  series.  We  saw  in  §2.5  (see  equations 
(2.45/  (2.46)  above)  that  equations  of  this  type  with  the  value  of  X 

prescribed  arise  in  the  analysis  of  certain  physical  problems.  In  what  follows 
we  shall  always  assume  that  the  value  of  the  constant  X  is  among  the  physical 
data  of  the  problem.  We  shall  follow  the  treatment  given  in  Srivastav  (1963  b). 


7.1 .  The  case  in  which  p  =  +1 . 

We  begin  by  considering  the  pair  of  equations 


OO 

^  n  an  cos(nx)  =  f(x),  0 

n*i 

00 

i  a„  +  ^  an  cos(nx)  =  °» 


<  x  <  c 


C  <  X  <  V  . 


(7.1) 


(7.2) 


Integrating  the  first  of  these  two  equations  with  respect  to  x  from  0  to  x  we 
obtain  the  relation 

^Xaox+2_/  an  sin(nx)  =  ^(x)>  0  <  x  <  c  (7.3) 

n*i 
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rx 

where  7kxj  =  '  fft'dv 

d  f* 

If  *»nsu',«  that  for  0  «•  x  <  r* 

•i# 

r*n 

a  ccr.(r  /)  = 


*  6  .  +  / 


#  /  '  - 

/ _ 

*  X  A.  (vf  ~  X,V> 


(7  4) 


thcr.  ^vpx.  f  mrulae  for  <j*+era> !  ning  the  Foxir.5  er  cceff  .lolenr-a  of  a  ha  If HTa.nfl!'* 
series •  ««■  hav*  the  equa+'  :ns 

a  fc  «  ,i  r%cMnryJ° 

l"*  _  :  f  2  S ‘  /  J  /  /  *  «?\ 


r  «  o 


j  X  ‘  U*  -  x’  T  i  O  'X  h/(u*  -  x 


Int.*r?hsr.«5vn«r  the  orders  of  the  Ir^egr&l-torw  *b*@*  <3  .-’/bis  integrals  we  cVvain 
the  ojrpreasx  ora 


a.  s 


g' vjdi;  . 


a.  ^  =  j  g(u]  J  -mj ').']d 


(’  5) 


If  re  sv‘*.w ti'.  j  .e  fbe^s  r.Kpr*»a.««ionn  an/!  ■>  equation.  (7  7 )  and  irvterobange  tbs 
order  of  tb»  auswatlon  sad  integration  «s  cbt-v  r  The  relation 


\1  S; 


^  n  'C 

g(u  *<  ,  .7  (nuj  nln(*'.x)fd'j  =  '?fx)  •  |  >  i ,  g(n)du,  0  <  t  <  o, 

J  o  -  /_*  ■■  *  j 

•*  • 

Haing  tbs  rpsu'j  ;.  f j  3<)  we  see  that  thin  *c\saticn  it  equivalent 

f  *  -JL<vteu_=  P(x)  ^  ,  T  t C  R(u)du  +  fC  g(v)du  r  Habfx^  x 

-c  //(*  -  u  .'•  ^  f>  *0  VO  sinhfffy) 

X  *  ,y  I,/uy)<5y  .  0  <  x  <  c., 

If  w  regard  -hi*  as  an  Integral  of  equation  of  Ateel  type  shoes  right-hand  side 
i*  a  known  function  cf  x  we  obtain  tbs  Fredholm  equation 


rc 

g{u) 

g(  ■  )  3?(  ■■  n/d  1  c  b/-.j) 

0  <  xi  ,  c , 

(7-7) 

where  fc(u) 

denotes 

tbs  free  term 

h(u) 

_  1  Jt 

?  -X?(x-  dx  __  ?  / 

(7,8) 

T 

«'  O  /•’  (U*  -.<■■')  r  V  o 

A/Vu*  -  X'; 

and  K(  t -  u)  denotes  the  kernel 


j-  oo  -  7T  y 

K(t,  u)  =  u  — ?--r-r - -  I  (ty)  I  (uy)dy-  ^  Xu 

'  -  '  J  smh  i.jryi  o  ?. 

•6 


(7,9) 


An  alternative  solution  leading  to  an  integral  equation  whose  solution  can. 
be  reduced  to  simple  quadratures  has  also  been  derived  by  Srivastav,.  We  put 


j  a n  +  /  a  cos  nx  =  co3'V2 


( £  x)  f  X  _ _ „  o  *  x  <  c 


J  c  a/(cg 


(7,1 0) 


CG3  X  -  COS  T. ) 


from  which  it  follows  that 


a,  =j—  j  ^t)dt- 


=  4^  J  <f>  (t)  j]?n(cos  t)  +  P^  ( co3  t)jdt,  (7.;1l) 


Substituting  these  expressions  in  equation  (7  -3 )  W5  obtain  the  relation 

OO 

—  X  f  (f)(t)dt  +  --  -  f  <1>(t)  }  jp(ccst)+P  (cos  t;j  sin(nx)  dt  =  P(x),  0  «x<c 

2\'2  Jo  '  V2  J  o  y}-  n 

from  which,  as  in  the  derivation  of  equation  (6.22)  we  deduce  that 


x  f  <f>(t)dt  +  cos Gg-x)  /*  ■  -  _  p(x),  0  ^  x  <  o, 

2/l/2  io  i o  \/(cos  t  -  cos  x) 

Inverting  this  equation  by  means  of  equations  (3  ,>4-2  a  and  b)  we  find  that 

2  d  f~  sin(-J  x)  P(x)  dx  1  d  r%  x  sin(-g-  x)  dx 

4>(t)  =  —  — -  I  - - —  Aa^  — /  “ - — 1  * -  -v 

7T  dt  /  0  /^(cos  X  -  COS  t)  7T  dt J c  ^/(coS  X  ",  cos  t) 

To  determine  s.  _  we  integrate  both  sides  of  this  equation  with  respect  to  t  from 
0  to  c  to  obtain  the  equation  ^ 

a  = - 1 -  (7,12) 


where 


1  +  X  i2 

V 2  j*  c  sin(^  x)  F(x)  dx  1  c  x  sin(-§-  x)  dx 

X  =--  -7 - ->  !,  =  -t  7 - -•  (7.13) 

7T  ,  0  \  COS  X  -  COS  C;  If  A/2  J  O  COS  X  -  COS  O ) 


Hence  we  have 
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2  d  - sin(i  x)  P(x)  dx  XI  d  x  sin(^  x)  dx 

(f){t)  =  —  —  I  V-  >  5-  1  !  ——————  , 

V  dt  Jt  V(  COS  X  -  cos  t)  Jr(l  +  XI  )  dt Ja  y(cosx-  cost) 


1,2-.  The  case  in  which  p  -  -1  „ 

We  now  consider  the  solution  of  the  pair  of  dual  series  relations 


2  X  a  + 
o 


oe 

rn  1  a  cos(nx)  =  f(x)*,  0  $  x  <  c 

J  " 

n-i 

oo 

O' 

ia-B  +  an  cos(nx)  =  0_,  c  <  x  $  n . 

n=i 

If  we  integrate  the  second  of  these  series  term  "by  term  we  find  that 

OO 

p 

^  a  ( ir  -  x)  -  )  n"1  a  sin  nx  =0,  c  <  x  <  w. 

0  ‘  LJ  n 

n-i 

We  now  assume  that  when  C  <  x  £  v 


oo 

^-X.aQ  +  V  n_1  a^  cos(nx)  =  sin(-J-  x)  f  .  -  .  - 

0  4—/  n  *  c  rJ (cos  t  - 

n*i  r 

which  is  equivalent  to  assuming  that 

a  =  “  f  [  f  (u)du  +  —  f  g(t)d  t  1 
T  Uo  z/2  Jc  J 

~  [  g(t)[p  (°°S  t)  -P  (cos  t)] 

A/2  J  c  Ln  n_1  J 

2 


1  d  t 

cos  x) 


an  =ntn  + 


dt 


where 


bn  =  ;  J 


f(u)cos(nu)  du. 


If  we  substitute  these  values  into  equation  (7,- 1 7)  and  interchange 
of  integration  and  summation  we  find  that 


OO 

J  g(t)S(x.,  t)dt  =  £  a s(ir  -  x)  -  l  b  n  s.in(n  x)  , 
c  n»=i 


c  <  x  i  v 


(7.14) 

(7.15) 

(7.16) 

(7.17) 

(7.18) 

(7-19) 

(7.20) 

(7.21) 

the  order 
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where 


S(x„.  t)  =  )  P  ( ccs  t)  -P  (cos  t)  sin(nx) 

a/2  U  I  n  n-’  J 


A/2U  Ln  n”  -J 

'  n  =  i 

sin(-^  x)  H(t  -  x) 

A/ (cos  x  -  cos  t) 

so  that  g(t)  is  the  solution  of  the  integral  equation 


sin(i  x)  I 


J x  a/(cos 


X  -  cos 


=  )  b  sin.  nx  -  a  (  w  -x) 


(7,22) 


with  solution 


1  d  r-e  r  a  d  e  ir  (  v  -  x)ccs(-g-x)d  x 

g(t)  = - \  b  { P  (cos  t)  +  P  (cos  t)|  - - j  - 

a/ 2  dt  L-i  n  *■  n  n  '  7T  dt  J  t  /l/(cos  t  -  cos  x) 


QO 

V'b^  |p  ^(co3  t)  +  P^  (cos  t)l  =~  [  f(u)  ff(up  t)  d  u 

C.—'  W  j  o 

nr. « 

•where 

oo 

<j(u;  t)  =  ^  ^P^  ( cos  t)  +  ^  (cos  t)'!  cos(nu)j,  (0«u<c<t<ir) 


f^2  cos(-g-  u)  H(t  -  u) 
Vi  cos  u  -  cos  t ) 


and  hence 


or  'J  2f/2  ,c  f(u)  cos(^u)du  2  rc 

)  ,bn  [Pr>os  +Pn-,(cos  *>]  =  -  /  - -  /  *00  du 

w  J  v  J  o  Vicos  u  -  cos  t)  wJo 


showing  that 


— - Y  b  fp  (cos  t)  +  P  (cos  t)l  =  -  —  [ 

^2  dt  L  "l"  n”  J  IT  dt  J, 


2d  *  c  f(u)  cos(^  u)  du 


o  //( cos  u  -  cos  t  ) 
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The  solution  of  the  integral  equation  (7.-., 22)  can  therefore  be  written  in  the  form 
2d  .c  f(u)  cos(-|-  u)  du  aQ  d  ^  (  ir  -  x)  cos(-g-  x)  dx 


a  f 

g(t)  = - [ 


IT  dt  Jo//  (cos  U  -  COS  t  )  7T  dt  it  >/{  cos  t  cos  x) 


C7.2J) 


The  value  of  a  can  be  found  by  a  simple  integration..- 
7., 3r  Equations  of  the  Second  Type,, 

If  we  make  the  substitutions 

c  =  tt  -  y  .  x  =  7T  -  y,  Xa  =  A  ,  X-1  =A5  n(-l)na  =  A  ,  f(n-y)=t  (y) 

0  0  n  n  1 

in  the  equations 


00 

4rXa  +  )  n  a  cos(nx)  =  0. 

2  0  L  n 


n  =  i 


0  ?  x  <  c 


OO 

£  a°  +  ^  a„  cos(nx)  =  f(x)..  0 


<1  <  i 


(7-24) 


(7,25) 


n»  1 


we  find  that  they  take  the  form 


OO 

£aAc  +  ^  n  " 1  A  ^  cos(ny)  =  f^(y).  0  «  y  <  y 


n*  1 


v  . 

0  n 

n  =  i 


^  A_  +  )  A  cos(ny)  =  0, 


y  <  y  <  v 


which  have  been  considered  in  §7.2  above  From  equations  (7-19) 3  (7„20)  and 
(7. 23)  we  can  write  down  the  solution  of  these  equations  and  reverting  to  the 
original  variables  we  find  that  the  solution  of  equations  (7, 24)  and  (7.25)  wn 
be  written  in  the  form 


Xao  =  W  [fa  fQ  g(t)dt+/c  f(x)dX] 


=  ~  j,  cos(nx)d x  +  -^~j  x 

|  P  (cos  t)  +  P  (cos  t)lg(t)d 
t  n  n-M  J 


(7.26) 


x 


where  g(t)  is  given  by  the  equation 


^  ^  2d  i’  it  f(u)  sin(4  u)  du  XaQ  d  px  sin(4  x)  dx 
v  dt  Jc  //(cos  t  -  cos  u)  7T  dtj04/(cos  X  -  cost) 


(7,27) 


In  a  similar  way  we  can  reduce  the  solution  of  the  pair  of  dual  series 
equations 


p  X  a  + 


00 

Zn”1  a  cos  nx  =  0, 

- 

n:-1 


0  <  x  <  c 


(7-28) 


r' 

4a  +  )  a  cos  nx  =  f ( x) c  <  x  s  ir 
0  /_/  " 


(7:29) 


n -i 


can  be  reduced  to  the  type  considered  in  §7-,1  above-  In  this  way  we  find  that 
the  solution  of  this  pair  of  equations  can  be  written  in  the  form 

.  7T 


a  =  [  g(t)  dt 

0  X  Jc 


(7.30) 


=  f  g(t)[p,(cos  t)  -  (cos  t)"|dt  (701 ) 


,-u  cos  (  4  *)  dx  d  .  *r  (  ir  -  x)  cos(4  x) 


2d  fir  r-U  coal  p-  XI  O.X  a  r  ■. 

g(t)  =  -  —  f(u)  du - a—  - -  dx. 

ir  dt-  Jt  J  t|/(cos  t  -  cosx)  °dt  Jt#/(cost-  cos  x) 


(7.32 


7; 4.  Cosine  Series  analogous  to  Fourier-Bessel  Series 
We  now  consider  dual  series  equations  of  the  type 


OO 

d 

n«  i 


(n  -  4)p  a  n  cos(n  -  4)  x  =  f  (x) . 

an  cos(n  -  4)x  =  fg  (x) 


I; 

n«i 


0  s  x  <  c 


c  <  X  ?  JT 


(7.33) 


(7*34) 


0  <  X  <  5 


(7,52) 
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ao 

£ 

n  =  i 


(n  ~  g)  a  coa(n  -  g)  x  =  f(x) o  <  x  «  f 


(7.53) 


can  be  written  in  the  form 


“  /  (1  +  cos  u)g(u)du,  a  =  -J—  /  g(uX  P  (cos  u)  -  P  (cos  ui  du,7 

V2  Jo  "  V2Jc  4"  3 


(n  >  2)  (7.54) 


where 


h/2-y  d  .-it  du  i  a  .  i 

g(t)  = - i  -  - — - 

dt  Jt  a/ (cos  t  -  cos  u)  ,  v  dt  Jt  //(cost-cosu 


Id  rn  du  -7T 

-I  f(x)dx  (7*55) 


with 


so  that 


fi 

/ — 

b 

1  fV  du  f1 

l  „  . 

1  c  //(cos  c  -  cos  Jl)  J 

vie  //(cos  c  -  cos  u)  Ju 

f(x) 


dx. 


For  this  solution 


g)  x  =  sin 


g(u)  du 


c  //(cos  u  -  cos  x) 


C  <  X  <  W. 


(7=56) 


The  study  of  dual  series  relations  in  which  the  series  involved  are 
expansions  in  terms  of  associated  Legendre  functions  have  been  studied  by  Collins 
(196i).-.  whose  treatment  we  shall  follow  here,. 

We  begin  by  considering  the  problem  of  the  determination  of  the  constants 


M 


such  that 


OO 

£% 

n*o 


T-  (cose)  =  f(e) 

m*n  '  ' 


0  *  6  <  a 


(8,1) 
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ao 

^  ( 2n  +  2m  +  l)anTm  +  n(  cos  9  )  =  g(  d)  u  a  <  9  «  ir 


(8-2) 


where  the  functions  f(g)  g(  g)  are  prescribed  and  we  assume  that  sin  ™  6  f(g) 
and  sin  m  9 gig)  are  continuous  in  the  closed  intervals  £o ,  aQ, 
respectively  In  §3.4  we  saw  how  dual  series  equations  of  this  type  can  arise 
in  the  analysis  of  problems  in  mathematical  physics  - 
We  make  the  assumption  that. 


i> 


+  2m  +  l)a  T  ^  (cos  q)  =  h(g).  „  (.0  $  g  sj  a) 

n  m+n 


(8.3) 


then  using  equations  (3  48)  and  (3  49)  we  find  that 


a  =  -5-' 
n  2 


=  i(-l)m  /  h(x)Tm  (cos  x)sin  x  dx  *  f  g(x)T|n"n(cosx)sinx  dx,;  (8-4) 

i  o  +  Jet 


(n  =  0V  1,  2  ■  .,;.).  Substituting  this  expression  for  an  into  the  series  on  the 
left  hand  3ide  of  equation  (8  l)  and  interchanging  the  order  of  summation  and 
integration  -  we  find  that  h(g)  satisfies  the  integral  equation 

,•  a  ra  r  v 

i(-l)  /  h(x)S  (g.  x)sin  x  dx  =  f(g)-i(-l)  /  g(x)S  (g  .  x)sin  dx  0  <  9  <  a  (8.- 5) 

Jo  m  Ja  m 


where  the  kernel  Sm(fl3  x)  is  defined  by  the  equation 

OQ_ 

^  m*n  (coa  T  m»n  (°°B  x)- 


(8.6) 


Now.  by  the  addition  theorem  for  Legendre  polynomials  (MacRobert;  1947» 

P  328),  we  have 

r 

Pr(cos  0)  =  Pr(cos  9  )Pr(cos  x)  +  2  (-1  )m  coa(m ift)  T^cos  9  )T™(cos  x)  ,  (8,7) 


where  r  is  a  positive  integer  and 


cos  0  =  cos  9  cos  x  +  sin  6  sin  x  cos  ip  , 
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From  this  it  follows  that 

\2  tt(-1  )  m  T"m  (oos  6)  T*(cos  x)  ,  r  *  m 

r2»  .  J 

j  P^(  cos  8  )  cos  m  iji  d  tfi  = J 

0  <  r  <  m 

and  hence ;  from  equation  (8  6 ) ,  that 

Sm(  e ,  r)  =  r 21 

Pm  +  n(cos  9)cosm$rdy  =  L±uy  r2  Pn  (co3  Q)cos  m  tjj  d p., 
1  2  7 r  l—1  Jo  *  2i  J  o 


Using  the  fact  that 


oo 

I 


P  (cos  8  )  =  (2  -  2  cos  9  )  2 


we  find  on  interchanging  the  order  of  integration  and  summation  that 
\m  r  2v  ___/_  ,n,,  /  .  \n  r  2v  — 


;  (  e  ,  x)  =  f  *  .ooaln  fe)  .  _  IdJL  f  U _ L 

m  2  v  J  o  tf(2  -  2  cos  8)  2v  J  o  V(s^  +  s2  - 


2s  s,  cos ' 
2  1  2 


where  s^  =  2  sin  g  6  cos  g  x;  s2  =  2  sin  g  x  cos  g  0  and  3^  05  s2  >  0  for  all 
8  and  x  since  both  lie  in  the  open  interval  (0-  n)..  The  integral  occurring 

in  the  equation  on  the  extreme  right  can  he  put  in  another  form  by  using  a  lemma 
due  to  Copson  (l947);  we  find  that 

[■  ,)  _  2(-l)’  /  “^V*^  d„ 

*(S.S,W°  “  sZ)(s*  -  s2)} 


S 


(8,8) 


1  2' 


If  we  substitute  the  form  (8  8)  for  3  (0,  x)  into  equation  (8,6)  and  notice  that 

m 

when  0  <  x.  s^  <  sz  and  when  0  >  x,  si  >  s2,  we  find  that  h($)  satisfies  the 
equation 


sin~m  0 
v 


”  /•©  rS  ZW,  fCt 

/  h(x)sin,’n'  x  dx  /  2  ^ - 8  1,  +/  h(x)sin1_mx  dx 

Jo  W[Is;  -  s2)(s2  -  S2)j  Jo 

s 

i  fl2m 

o  //Qs‘  -  sa)(s*  -  S*)J 


L 


I 

1 

I 

I 

I 

I 
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,  -m  .  rs  2m 

=  g(x)aln’-mxdx  1  ^ - -r,  (0  <0<a).  (8.9) 

v  Ja  Jo  //fa2  -  a2)(s^  -  s2)j 

Ve  now  change  the  variable  in  the  inner  integral  from  s  to  u  where 

s  =  2  cos  g  6  cos  g  x  tan  g  u 

and  invert  the  orders  of  integration  in  the  integrals  in  equation  (8„9)  to  obtain 
the  integral  equation 

pQ  H(u)  (tan  g  u)m  du  .©G^u,  a  )  (tan  g  u)2mdu 

/  - =  (tan£fl)‘  f(0)  -  I  — - - — —  (0  <  0  <  a)  (o01 0) 

Jo  V  ( COS  U  -  COS  s)  Jo 


V( 


COS  u  -  cos 


e  ) 


where  the  functions  H(u)  and  G(u,  a  )  are  defined  by  the  equations 

1  i* a  h(x)  (cot  g  x)"  sin  x  dx 

H(u)  =  —  / 

2  v  Ju 


A/( 


coa  u  -  cos  x 


) 


1  rir  g(x)  (cot  g  x)m  sin  x  dx 

G(u,  a)  = —  /  — - - -  (u  <  a). 

2  v  Ja 


V( 


cos  u  -  COS  X. 


) 


(8.11) 

(8.12) 


Prom  equations  (3  ,42  a  and  b)  we  see  that  the  solution  of  the  integral 
equation  (8,1 0)  is 

(cot  g-  u)m  d  i‘ u  sin  0 ( tan  g  0  )m  f(fl)de 

H(u)  =  —————  —  / - G(u,  a  )  „  (0  <  u  <  a)  (8*13) 

it  du  J  o  cos  9  -  cos  u) 

and  from  equations  (3,43  s  and  b)  that  the  formula  giving  h(x)  in  terms  of  the 
function  H(u)  so  determined  is 


h(x)  =  -  2(tan  i  x)m  cosec  x  ±  f “  sin  u  du  ^  <x<a)o 

dxJ  x  V(  cos  x  -  cos  u) 


(8.14) 


The  corresponding  value  of  a^  is  obtained  by  inserting  this  form  for  h(x) 
into  equation  (8,4)0  However,  it  is  possible  to  determine  a  in  terms  of  H(u), 

Collins  has  shewn  that 

a  _  (2i)  T  v,n  +  2ra  +  1  )f  f  (ain  £  u)  cos(g  u)  II  (cos  u)du 

r(n  ♦  m  +  %)  Uo  m>n 

+J  G-qC'i)  (sin  ?  u)"  cos(g  u)nm  n  (cos  u)duj 


(8.15) 
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is 


if  sin  n  a  sin(n  +  l)a_  ^  (8,25) 

n  7T  L  n  n  +  1  J 

The  corresponding  expression  for  H(u)  is  *-  cos  ^  u  so  that  when  0  i  9  <  a 
we  have  from  equation  (8  <,22)  that 


i<2n 


+  1  )  a  P  (cos  6  )  = - cosec  6  — 


2/^2  &  ra  sin  u  cos(-g;  u)  du 


tt  ddJS  (cos  9  -  cos  u) 

n*o 

The  integrations  are  elementary  and  we  find  that 


oo 


+ 1 )a  P  (cos  s) 


ir 

=  —  j  ^  7r  -  sin"1 
7 r  ^ 


/cos  -5-  a  \  cos  -j  a 


icos  ifil  ^/(cos^  ^0  -  cos2^a)- 


,0  s6<  a.(8„26) 


Similarly  if  we  put  m  =  i  in  the  general  formulae  we  find  that  the  pair  of 
equations 

OO 

Za  T  ( cos  6  )  =  f(s)  ,  0  $  6  <  a  , 

"  n  +  1 


(8.27) 


n  =  o 


oo 

^  (2n  +3)  a  n  Tn'(  (cos  8)  =  g(fi);  a  <  6  s  w  , 


(8,28) 


n*o 


possesses  the  solution 

a  =  -1— e  f  H(u)  tan(4-  u)  sec(|-  u)  f(n  +  1 )  sin(n  ♦  2)u  +  (n  +  2)sin(n+l)u]  du 

n  V£Jo  L  J 

/•  7T  ^ 

+  J  G-o(u)  tan(-g-  u)  sec(^  u)|(n  +  l)  sin(n  +  2)u  +  (n  2)sin(n  +  l)ujdu  (8,29) 


where 

2  cot  75-  u  d  r  u  sin  (■§■  d)  f(fi)  d0 
H(u)  =  — — ——  —  j  1  >  *■  G-(up  a  ) 


1  r* 

G(us  « )  =  “ 

ir  Ja 


m  du  J  o  /^(cos  0  -  cos  u) 
1  r*  s(0  COS  *  2  (i  e)  de 


V (cos  8  -  cos  u  ) 
and  G-  (u)  is  defined  by  equations  (8  =,17) 


(8.30) 


(8,31) 
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We  now  consider  the  equations 


00 

E 


n  n  m-;  n 


(cos  6  ) 

ii 

Hi 

•25 

V 

(0  <  Q  <  a  )  ; 

(8.32) 

(cos  6  ) 

=  g(e)_ 

(a  <  6  <  v  ) ; 

(8,33) 

where  H  n  is  a  given  function  of  n  which  is  0(n  ‘)  for  large  values  of  n.  As 
before  we  assume  that  sin~mef(e)  is  continuous  in  (Or  a  )  and  that  sin^  6  g(o) 
is  conti minus  in  (a  s  tt  ) „  If  we  make  the  assumption  (8r3)  then  the  coefficients 
an  are  again  given  by  equation  (8... 4);  if  we  insert  the  value  (8,4)  for  a  into 
the  series  on  the  left  hand  side  of  equation  (8„32)  we  obtain  the  equation 

OC^ 

\  T (cos  e)(l  +  Hn)  f  h(x)  T™in  (cos  x) 


sin  x  dx 


=  f(e)  -^(-l)m^  (l  +  H  )  T-m  (cos  e)[  g(x)Tra  (cosx)sin  xd  xn(o<  fi  <aX(8c34) 

/  t  n  m+n  j  ^  ni4-n 


n*c 


As  in  the  derivation  of  equation  (8„ 1 0)  we  find  that 

•ff  r®G(u;a)tanZ,I^u  du 


5-(-l)  T”^  (cos  d)  [  g(x)Tm  (cos  x)sin  x  dx  =cot rg-0  / 

Li  m>n  Ja  m:n  Jr- 


^-(-if  T  X  (cosd)f  h(x)Tm  (cosx)sinxdx  =cot / 

L_/  m*n  J0  n'  +  n  J0 


o  4/(cosu  -  cos  6  ) 
^H(u)  tanzm-g-udu 


^(cosu  -  cos  6  ) 


(0,35) 


(8.36) 


where  H(u)  and  G(uP  a  )  are  defined  by  equations  (8.1 1  )  and  (8,12)  respectively. 
If  we  write 

t a 

I  =  ^-(-l)  7  H  T  m  (cos  e  )  /  h(x)  T m  (cos  x)sinx  dx 

/  ,  n  m+n  j  in  +  n 

L—J 
n*o 

then,  using  equations  (3,47)  and  (3 -50),  we  find  that 


’T"'  H  r  (n  +  2m  +  1 )  r< 

=  i)  - — ; - —cotmie 

r(n+0  J  c 


®Rm  (u)  du 

m^n  '  ' 

cot 

)  Jo4(cosu-cos0) 


/: 


h(x)sinx  cotm-g-x  dx  x 

rx  R  ( v)  d  v 

/  •«, — sain. 

/  -  If 


4  o  4/(cos  v  -  cos  x) 


t 
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It  can  be  shown  that 


I  =  cot 


'£»  f 

Jo 


9  tan  m  i  u  d  u 


A/(  cos  u  -  cos  6  ) 


Ja  tanra£  vK^u,  v)H(v)dv  (8.37) 


where 


OO 

v)  =  ircotffl^u  cotm^  Hn  R^n(u)  Rm"n(v).  (8.38) 


Similarly  if  we  write 

OO  Jf 

j  =  i(-l)m^  H  nT*"n  (cos  e)j"  g(x)  Tmm+n(c°s  x)sin  x  dx 

n  =  c 

it  can  be  shewn  that  if 


(8.39) 


&  (v)  = 


'G(v„  a  ),  o  <  v  <  a 
&(vv  v)  ,  a  <  v  <  ir 


then 


-6  tan  j  u  du  fir 

J  =  cotra  -  j  tan^vK^u,  v)Go(v)dv.„ 

-/  O  V(  COS  u  -  cos  e  )  -  o 

If  we  substitute  from  equations  (8„35)3  (8,36)_-  (8.37)>  (8, .41 )  into 
equation  (8.;34)  we  find  after  some  manipulation  that  the  function 


(8,40) 


(8„41 ) 


J(u)  =  tanm  ^-u^H(u)  +  G(u,  a  )J 


(8,42) 


satisfies  the  Fredholm  integral  equation  of  the  second  kind 

r  a  cotm -^u  d  r u  f(x)  tan^x  sinxdx 

J(u)  +  K  (u,  v)  J(v)  d  v  = - /  — - - - - - 

Jo  1  w  du  J  o  v(cos  x  “  003  u ) 


r  v 

-  /  tanra  j  vK  (u,  v)  Gr(v_  v)dvs  (0  <  u  <  a)  (8<-,43) 
J  a  1 


with  symmetric  kernel  K  ^  (  u ,  v)  defined  by  equation  (8,38),  Once  we  have  found 
the  solution  of  this  equation  we  can  determine  H(u)  from  equation  (8,42)  and  the 
function  h(x)  from  equation  (8,14). 


16 


The  case  m  =  0  is  of  special  interest..  We  then  have 

J(u)  =  H(u)  +  Gr( u  .  a  ) 

where  J(u)  is  an  even  function  of  u  satisfying  the  integral  equation 

J(u)  +'-\a  Ko(u  -  v)  J(v)  d  v  =  i  ±  r 

tr  .  ..a  Jr  du  Jo  a/(co&  x  -  cos  u) 

\  H  cos(n  +^)uf  G(v,  v)cos(n  +|)vdv.  |u.|  <  a 
ir  Lj  n  J  a  ' 

n^o 

where  « 

Kq(x)  =  ^  cos(n  +  x 


(8,44) 


(8,45) 


(S.46) 


9  Dual  Equations  involving  Series  of  Jacohi  Polynomials 

Recently  Noble  (l 9&3)  and  Srivastav  (l9^3^)  have  considered  dual  series 
relations  involving  Jacobi  polynomials.  In  this  section  we  shall  follow  the 
method  of  solution  due  tc  Srivastav, 

We  consider  the  equations 


p  T  ( a  +  n  +  £)  {a  o) 

)  — ; - "AnPn  (cos  e )  =  f(e);,  o  $  e  <  $ 

L~  r  (  a  +  n  +  1 ) 


(9,1) 


P  r(^  +n+;)  _(a , fi)  ,  , 

/  — : - : - - - -  A  P  (cos  8)  =  g(e).  6  $  TT 

yr  Tip  +  n  +1)  "  " 


(9=2) 


1  \*\P) 

vhere  a>  ~  %r.  P  >  ~  ^  and  P^  is  the  Jacobi  polynomial  defined  by  the  equation 


3(  a  ) 


(x)  =^n  *a^j  2?  (-nv  n  +  a  +  j8  +  1;  a  +  1  f  £  -  £x) 


(9=3) 


035  alternatively,  by  the  equation 


(*)  js  *  1*  4  »**). 


(9.4) 
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We  shall  also  require  the  fact  that  if  a  >  -1 ,  /3  >  -1  . 

(l  -  x)  a  (l  +  x)^  ^  ^  (x)P^“s  ^  (x)  d  x  =  a  '  ^  ^  6  , 


m  n 


(9.5) 


where 


(  §  )  _  2a  +  ^  +  ^  ?  Ca  +  n  +  1 )  P  (ft  tntl) 

’  n  nj  T(a  +  /3  +  2n+l)r(a+ft  +  n  +  1 )  : 


(9o6) 


We  observe  that  if  the  solution  of  the  pair  of  dual  equations  (9-1 )  and  (9.2) 
is  known  we  can  derive  the  solution  of  the  pair 


OO 

I 

n»e 


r(  a  +  n  +— )  / 

- 8 —  a  (cos  e)  =  f(e) ,  0$  e  <  4> 

r(a  +  n  +  i)  n  n 


(9*7) 


co 

I 

nrttf 


r(/9  +  n+i)  / 

- a  p^  ,/3'  (cos  e)  =  g(e)„  <t>  <  e  $  n 

r  (/9  +  n  +  l)  n  n 


(9.8) 


by  a  simple  change  of  variables  >  Por  that  reason  we  shall  consider  only  the  pair 
(9,1)  and  (9.2). 

If  in  equations  (8,.l)  and  (8.-2)  we  make  the  substitutions 

T  ”m  (cose)  =  X„(P  +  ijsinli.  p(m  ,  m )  (cQS  Q  ) 

2  m  r(n  +  m  +  1)  " 


r  (m  +  n  +  £) 

-  A  ,  f(e)  =  2  sin  e  P(e) 

r(n  +  l)  n 


g(fl)  =  2  1  "m  sin m  8  die) 

we  find  that  these  equations  can  be  written  in  the  form 


OO 

l 


OO 


n-,o 


r(m  +  n  +  ■§•) 
r(m  +  n  +  1 ) 

r(m  +  n  +  ^) 
T  (m  +  n  +  1 ) 


A  P^ 
n  n 


m  ) 


^  p(  m  p  m  ) 
n  n 


(cos  e )  =  f( e)0 


(cos  e )  =  G-(e) 


0  s  6  <  a 


a  <  6  $  it 
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so  that  the  equations  (8.1 )  and  (8„2)  are  special  cases  of  the  equations  (9.1)  and 
(9.2)  with  a  =  /3  =  m„ 

We  begin  by  considering  the  integral 


u  sin  x  P  a,J  ^  (cos  x)  (sin  £  x) 


2  a 


/✓(  COS  X  -  cos  u) 


■dx3  0  «  u  5  ir  ,  (9.9) 


If  we  replace  the  Jacobi  polynomial  by  the  hype rgeome trie  function  given  by  equat¬ 
ion  (9-3)  and  interchange  the  order  of  summation  and  integration  we  find  that 


I  (u)  =  2”  “tt2  +Jl  _  cos  u)“+2'Pa+  2;  P-  2  (cos  u)#  (9,10) 

r (a  +  n  +-) 

Similarly  we  can  show  that  the  integral 


!,(▼)  = 


n  sin  x  P  @  (cos  x)  (cos  x)^^  dx 


V( 


COS  V  -  COS  X, 


) 


(9.11) 


has  the  value 


(9,12) 


I  (v)  =  2“%  2  (l  +  cosv/+2  r  +.n,+.ll  p  «-  +  2  (cos  v)  . 

2  r(/?+n+-) 

Prom  equations  (9o10)  and  (9,1 2)  it  follows  immediately  as  a  result  of  simple 
differentiations  that 

—  I  (u)  =  2~a  7 r  2  ’*  1  +  ^  (l  -  cos  u)a~  2  sin  u  Pa  2  >  P  +  2  (Cosu) ,  ( 9 1 1 3 ) 

du1  r(a+-g-+n)  n 

and  that 

-  —  I  (v)  =  2~  ^  w  2  ^  ^  +  -■■+■^1  (1  +  cos  vf  2  sinvPa+^?^  2  (cosv).  (9&14) 

dv2  r(^9  +  i+n)  n 

Using  the  orthogonality  condition  (9o5)  we  can  easily  show  that  if 

*  & 


f  (cos  v)  =  ^  cn  p“  '  (cos  v)  ,  0 


a  V  i  V 


(9.15) 


n*o 


then  the  coefficients  c  are  given  by  the  equation 
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.  =  Sal (p  +  0+  2n  +  l)r(g  +  jS  +  n-t-  1 )  f "  (ainjv)2a +1  (coa£ 
n  r(a  +  n  +  1  )r  (/9  +  n  +  1 )  Jo 


•g-v)^  f(cos  v)  x 


x  P^n  *  ^(cos  v)  d  v  .  (9-1 6) 


We  are  now  in  a  position  to  solve  the  pair  of  equations  (9d),  (9.2).  If 


we  multiply  both  sides  of  (9-1 )  "by  sin©  (sin  g©)^a(cos  ©  -  cos  u)  2  and. 
.integrate  with  respect  to  ©  from  0  to  us  then  differentiate  with  respect  to  u, 
we  find  on  using  equation  (9.13)  that 


00 

i  u)a  2  sin  u£  A^pJ*  2’/3+2(cos  u)  =  P(u)_  0  <  u  ^  <p 

n-r, 

where  P(u)  is  defined  in  terms  of  f (© )  by  the  equation 

P(u)  =  —  f  sinfi  (sin  g©  )^a  (cos©  -  cos  f(©)d©„ 

J  o 

Similarly  if  we  multiply  both  sides  of  equation  ( 9-2 )  by  3in©(cos  g©)^ 
(cos  u  -  cos  e)-2  and  integrate  from  u  to  w  we  find  that 


2_r"  jr  2  (l  +  cos  u^f  A^P^3  2j^  +  2,(cos  u)  =  G-(u)  , 

r.  =  o 


(9.17) 


(9.18) 


(9.19) 


where  G-(u)  is  defined  in  terms  of  g(©)  by  the  equation 

Gr(u)  =  j  sin  ©(cos  g©  ffi  (cos  u  -  cos  0  )"^g(0)ae„  <t>  <  u  s  j r 
J  u 

If  we  now  make  use  of  the  formula  (9.1 6)  we  find  that 

A  =  (2g  )~2  nl  (  g  +  0  +  2n  +  j)  V  (  a  +  jjULjtJil  x 
n  r(a  +  n+g)r(/S  +  n+^) 

x  j(cos  g$)^ 1  P3  2  9  ^  +  2  (cos  <f>)J  sin  ©(sing©)^a(cos©  -  cos<£  )  2f(©)d© 

/•<#>  1 

+  (n  +  /?  +  g)j  sin0  (sin£  ©)  af(©)d©J ^  (cos  ©  -  cos  v)”2  (cosgv^sin gv  x 


(9,20) 


)a  +2*  &  ~  2 


2  (cos  v)  d  v 


8:0 


P 


a 


n 


i_ 

2 


(9,21) 


If  the  forms  of  the  functions  f  and  g  are  complicated  the  determination  of 
the  coefficients  Ar  by  means  of  equation  ( 9  ■.  21 )  is  obviously  a  very  difficult 
procedure 3  but  in  many  physical  problems  it  is  sufficient  to  know  the  behaviour 
of  the  functions 


ao 

f  (e)  =y 

P ( a  +  n  +  2)  h  ,  8)  ,  . 

(9,22) 

1  Lj 

n  =  0 

r(a  +  n  +  l)  n  n 

00 

8,(0  =  > 

r(/3  +  n  +4  )  (a./S)  .  \ 

■  - —  A  PN  (cos  0  )  *  0  <  0  <  <£ 

(9. 23) 

1  Lj 

n  =  o 

r(p  +  n  +  i)  "  " 

We  can  write 

equation  (9-17)  in  the  form 

y*  A  ^  P^1  2  ;  ^  +  'S'J^cos  u)  =  Za  v  2  (l  -  cos  u)2  acosecuP(u)# 

n  =  o 

If  we  multiply  both  sides  of  this  equation  by  (sin  u)(sin  u)2a  1  x 
(cos  u  -  cos  x)~2  .  4>  <  x  £  it  and  integrate  we  obtain  the  relation 

2i~avi  (i  -  cob  ef  f  (e)  =  £  f  — Fi'ilAv - 

Jo  tf(  COS  U  -  COS  6  ) 

Q 

+  (2b-)  sin  u(sin^-u)2<1  1  (cos  -g-  u)~2^  "*  - G;(u)  du 

J  $  V(cosfi  -  cos  u) 


i,er.  when  <f>  <  8  «  v 


f  (e)  =  2  a  b  (l  -  cos  6)’a  [  - pfel-du 

1  ■'o  //(cos  u  -  cos  8  ) 


+  2a  tt  ^  (l  -  cos  8 )  af  sin  u(sin-g-u)2a  (cos-g-  u)~2^  — —  ,  (9.24) 

J  4>  ^(pozd-ccB-u) 

where  F(u)  and  G-(u)  are  defined  by  equations  (9d8)  and  (9?20)  respectively. 

Similarly  f  or  0  «  8  <.  <f>  we  have 


f  (0)  =  -  2P  v  ^ (  1  +  cos0)  ^  cosec 

2  &eJe 


d  r  4>  (cos-g-  u)^sinu(sin  g-u)~2a 


d 6J8  //(cos  8  -  cos  u) 


P(u)du 
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.j 3  -1/..  . \-fl  .  d  fV  sin  u  G(u)  du 

~2r  7T  (1  +  cos  e)  'oosec 0— /  - *  — 

dai<£  f/(coad  -  cos  u) 


(9.25)  ' 


W. 


J. 


E. 


L. 


H. 


T. 


B. 


E. 


I. 


I. 
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G. 
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